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FOREWORD 


This report combines material generated under Contract NAS z~ 
8738, Aircraft Differential Maneuvering , during the period February 197t> 
through January 197(5. It consists of updated versions of three reports 
issued during the period, listed as follows. The technical monitors were 
Dr. Mark Ardsma and Mr. Michael Tauber, both of NAS A -Ames Aero- 
nautical Systems Division. 

AMA Report No. 75-27, CALCULATION OF DIFF ERENTIAI .-TURNING 
BARRIER S UR FACES 


AMA Report No. 75-26, A USER’S GUIDE TO THE AIRCRAFT ENERGY* 
TURN AND TANDEM-MOTK 


AMA Report No. 75-7, A USER'S GUIDE TO THE AIRCRAFT ENERGY- 
TURN HODOGRAPH PROGRAM 
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SUMMARY 


The computation of composite differential-turn trajectory pairs is studied 
for "fast-evader" and "neutral-evader" idealizations introduced in earlier publica- 
tions. Transversality and generalised corner conditions are examined and the 
joining of trajectory segments discussed. A criterion is given for the screening 
of "tandem-motion" trajectory segments. Main focus is upon the computation of 
barrier surfaces. Fortunately, from a computational viewpoint, the trajectory 
pairs defining these surfaces need not be calculated completely, the final subare 
of multiple -subarc pairs not being required. Some calculations for pairs of 
example aircraft are presented. 


INTRODUCTION 


The modelling of a turning cha.se in air combat as a differential game with 
on ergy-mo defied vehicles is described in Ref. 1 and the characteristics of families 
of solutions are investigated in Ref. 2. The present report, a sequel, deals with 
the computational solutions of families, with emphasis upon the determination of 
barrier surfaces, those surfaces that separate successful pursuit from successful 
evasion. It will be assumed, for brevity, to be read in conjunction with Refs. 1 
and 2. 
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TRANSVKRSALITY CONDITIONS: ' ’ FA ST -EVADE El” MODEL 


A simplified vers Lon of the problem is of interest initially for illustrative 
purposes. Thus attention is directed to the •’fast-evader” altitude-dynamics ideali- 
zation, in which both pursuer and evader have a free eholce of altitude within their 
respective bounds and no intervals of "tandem motion”, with the pursuer driving 
the evader upward, appear (Ref. 2). 


Consideration of a related problem with vehicle models further reduced in 
order (Ref, 2) suggests that fast -evader trajectory pairs have the following form. 
For the aircraft having the higher maximum sustainable „ara rate as evader , there 
is simple pairing of Euler solutions terminating in the "capture, set", a subset of 
the "target set" (both defined in Ref. 1) or on its boundary. If the sustainable-tum- 
x-ate -superior aircraft is pursuer, and if it is inferior to the evader in maximum 
instantaneous turn rate, some of the trajectory pairs have two subarcs. The initial 
subarc almost attains turn-angle closure with loft-ceiiing-matcked energies; the 
second comprises a high-turn-rate pair of trajectories joined to the first at a cor- 
ner and terminates in the capture set or on its boundary. 


If the capture set lies entirely above the 1 o f t- e e i 1 mg ™ mat c h curve in energy 
space (pursuer specific energy chosen as the ordinate), capture necessarily takes 
place with the loft-eeiling inequality met with a margin, and both energy multiplier 
variables vanish at the final point by trans vers alitv. Such a situation arises in the 
case just mentioned exhibiting evader superiority in maximum instantaneous turn 
rate at matched loft-ceilings, as in Fig. 7 of Ref. 1, which applies for aircraft B 
pursuing A, 
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If the loft-ceiling-match curve forms part of the boundary of the capture 
set, as In Fig. 6 of Ref. 1 { A pursuing B ), the energy multiplier final values 
for capture along the loft-ceiling-matcli portion will be related by transversality. 

With h_ and h tile loft-ceilings for evader and pursuer, respectively, the 

L 1 2 
inequality 


V - V 2 0 


< 1 ) 


.and the equation 0 ? together defining capture, may be treated formally by 

adjoining their left members to the criterion function, final time t^, to form the 
augmented function: 


3? t f + A A Xf + A h < h L " h L > 

i-f 

f f 


(2) 


The usual transversality analysis leads to 
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(3) 

(4) 

(5) 


Imposing these conditions at a selected terminal capture point on the loft-ceilings 

match curve in energy space defines a one-parameter family of trajectory. pairs, 

i.e, , one of the two parameters A ^ and A. is defined by (3), while the other 

X 1l l 

can be regarded as the parameter of the family. 
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The limiting trajectory pairs corresponding to near-miss or grazing trajec - 


tories are of special interest and correspond to H = 0 . Such pairs are not proper 

members of the family but define the barrier surfa ce separating successful evasion 

and successful pursuit. The special values of the ratio A /A, that make II 

X h L 

vanish may be determined iteratively. There may be one, two, or more of such 
values, or none. Some insight into this is provided by the relationship between 
the dual hodograph figures introduced in Ref, 1. Such a figure is drawn for loft- 
ceiling -matched energies using y for the ordinate and {dh^/dE) E for the 
abscissa. The number of intersections of the two hodograph figures in the upper 


quadrants equals the number of vanishings of II . If one figure lies entirely within 
the other, there are none; this means one aircraft is locally superior and that there 
are no grazing trajectories through the point. An intersection in the upper quad- 


rants implies 11 = 0 by the existence of a common tangent to the figures; a per- 
pendicular to this through the origin defines the direction of the corresponding 
multiplier vector. The term "tangent” is employed loosely here to include contact 
at a corner of the figure. 
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COMPOSITE "FAST-EVADER" TRAJECTORIES 


If the pursuer does not have any land of low-energy superiority, the trajectory 
pairs are unbroken {except for what might be called ’’ordinary” corners, those satis- 
fying the Erdmann condition) and can be found by integration of Euler solutions back- 


wards in time using final multiplier values determined from eqs. (3)-(G); trajectories 


terminating at various points of the capture set's boundary computed in this way deter- 


mine the barrier surface. 


What might be called generalized corners, i. e. , corners with multiplier jumps, 
appear when the pursuer is superior in maximum sustainable turn rate and at the same 
time inferior in maximum instantaneous turn rate along the matehed-loft-ceiling curve 
in energy space. When the opposing' aircraft have this relationship, the pursuer, even 
having attained near angular closure with sufficient energy, is unable to effect capture 
immediately but only manages to force the evader into a high-turn -rate maneuver mid 
to steer the chase toward lower energies where capture may take place. It can be 
argued that sustainable-turn-rate-superiority of the pursuer along the matched-loft- 
eeiling curve to the lowest energies permitting flight makes closure with matched-loft^ 
ceiling energy tantamount to subsequent capture. The appearance of at, least some such 
two-subarc trajectories is insured with less overwhelming superiority, L e. , superiority 
of pursuer in maximum sustainable turn rate with energies of both craft unrestricted; 
in this case, a segment of the matehed-loft-ceiling curve will be the locus of generalized 
corner junctions. The upper end of this segment is defined by the limit of the pursuer's 
sustainable-turn-rate superiority. The lower end will be on or below the lower edge 
of the band of pursuer energies defined by pursuer sustainable turn rate superiority 
over the evader's highest sustainable turn rate detei-mined with his energy unrestricted. 
The arguments needed to establish that closure with matelied-loft -ceiling energy along 
this segment insures subsequent capture are elementary and similar to those employed 
for the sufficient condition of Ref. 1. 

Attention is now directed to the generalized corner junction of extremals appearing 
in this case where the pursuer has superiority at low energy but lacks a margin of maxi- 
mum instantaneous turn-rate to turn closure immediately into capture at high energy. 


Consider a family of Euler solutions through an Initial point having 0 and pur- 

suer energy high enough that some of the family members intersect Ax = 0 above 
the lost-ceiling-match curve. Call T(E ,.E ) the time -to -capture measured from 

X 

arrival at A\ = 0 , and At (not necessarily small) the time elapsed during the 
motion along the first subarc. Then an augmented performance index may be set up as 

P - T(I i )+*At+v(i -h T ) + A ' Ax (?) 

12 1-, X 

where the superscribed bars denote values at arrival at the surface Ax~0* 

If the inequality 



is met with a margin, then (i= 0 . The equality will obtain in many eases, however, 
since the pursuer will gent. rally tend to sacrifice energy margin for turn-rate in 
effecting closure. 


A transversallty analysis carried out with the inequality (8) and Ax = 0 
regarded as an intermediate target set leads to the system: 


H 

-1 

(9) 
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X 
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\ 

(10) 
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V 

5T h l 

- V dE 

1 

(11) 


dh 


\ = 

5T . L 2 

— -I- V 1 

hi dE 2 

(12) 


bE 2 2 





For extremal subarcs entering the surface Ax = G with strict inequality (>) in 
(8), v-0 and tlie conditions (9) -(12) determine all of the multipliers. For vj °> 
there is a one-parameter family whose limiting members occur for v unbounded; 
these are not proper members, as II - 0 . The ratio of multipliers 



is well defined in the limit. These multipliers and X _ are defined only within 

X 

a multiplicative constant for these grazing trajectories. (Note that this corresponds 
tc the classical definition of abnormality in variational problems: nonuniqueness of 
the multipliers.) 


The generalized corner condition of Bernhard (Kef. 3) furnishes a guide to 


which components of the multiplier vector jump and which do not, viz., the jump 
is normal to the singular surface. For i/«0* this surface is the portion of the 
plane Ax “ 0 above the matched-loft-c eiling curve, and only the multiplier X 

X 

jumps. For v 0 , it is the boundary defined by the equality sign in (8), and 


the energy multipliers jump to tlieir values bT/bE^ and bT/bE^ from the 
values (11) and (12); there is no jump in the component of the energy-multiplie r 
vector along the loft-ceiling-match curve’ s tangent. 


Constructing nongrazing trajectory pairs that have two subarcs and ter- 
minate in capture from the first-order necessary conditions is fairly laborious. 
First, "open-loop” trajectory pairs connecting the capture set and the inter- 
: mediate target set defined by (S) and Ax ~ $ are obtained. These are the second 
subarcs of composites and their energy multipliers provide the partials bT/ bE^ 

■ and bT/bE needed in (11) and (12). The system (9)— (12) then provides multi- 

i. .... ; cj ... ... 

plier end conditions for backward integration of Eul.er solutions which are the 
appropriately matched first subarcs. 



A stroke of good fortune is that grazing trajectory pairs, of main interest 
because they define the barrier, are relatively easy to calculate. This limiting 
case has the first-subarc terminal values independent of the second-suharc solution; 
the partials &T/6U and &T/bE do not appear in (13), which greatly simplifies 
the calculation of the barrier surfaces, A similar situation occurs with the neutral- 
evader model, to be taken up next. 


S 



' 'NEUTRAL- EVAD ER 1 r MODEL 


Adoption o-f the more realistic "neutral-evader" model lead s- to families 
of solutions having a more complex structure, as discussed in Ref. 2 . Under 
circumstances somewhat similar to those in which two-suoave solutions arise 
with the simpler model, solutions having a third, middle, subarc of so-called 
"tandem motion" appear. The tandem motion may itself be complex, including 
the singular arcs discussed in Ref. 2 . The discussion of tins section will ignore 
the possible appearance of singular sub arcs, whose role is not fully understood 
at present. It Is worth noting, however, that all of the singular subarcs found 
computationally in the examples to date have been screened out by one test or 
another, hence the assumed non-appearance may not be unduly restrictive. 


Trajectories of "tandem-loft" type, calculated by the procedure and com- 
puter program described in Ref. 4, are shown in energy space in Fig, 3 for air- 
craft B as pursuer and A as evader. There is one trajectory of particular 
interest that bifurcates, dividing the family. Those above it generally produce a 
loss of pursuer energy and originate at high pursuer energies. Those below the 
bifurcating trajectory correspond generally to increasing pursuer energy (although 
not monotonieally) , and originate at low pursuer energies. The two branches of 
the bifurcating trajectory form a separatrix , a boundary outside of which the. tandem- 
loft motion does not attain I oft -ceiling match. 

The family of tandem/loft trajectories' is determined by numerical integra- 
tion using a computer program such as that described in Ref. 4. As a preliminary, 
loft-ceiling rates at matched-loft-ceLJ ing altitudes may be compared at various 
points along the loft -ceiling -match curve in energy space, as in Fig. 5 cl Ref. 2, 

This was carried out with the computer program of Ref. 5, which also performs 
a single step of integration backwards in time using a coarse integration technique 
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(Euler extrapolation) to assist the search for the bifurcating trajectory and the 
separatrix. 


The region in energy space between the separatrix and the loft-ceiling- 
match curve plays a role in the neutral-evader trajectory family analogous to that 
of the loft-ceiling-match curve in that of the simpler fast-evader model. Trajectory 
pairs corresponding to pursuer closure with insufficient energy undergo transitions 
in this region to tandem -loft motion, which may continue until the loft -ceiling -match 
curve is reached: then, either capture occurs or a trails it ion to a high-turn- rate 
open-loop-optimal pair. The grazing trajectories for H = 0 join with tandem-loft 
arcs on the separatrix when the hoclograph figures are in the proper relationship, 
previously described; however, the scaling factors on the energy-rate axes of these 
figures are different: they are the magnitudes of the energy multipliers. The ratio 
of these is of main importance, and, along the separatrix, this ratio for H - 0 is 
obtained as 



This may be viewed as the determination of the normal to a two-dimensional trajec- 
tory from the components of velocity. The relationship (14) is employed like (13) in 
the determination of multipliers from t r an s vc r s ality conditions generally similar to 
those given for the fast -evader model. The computation of barrier candidates comes 
down to the use of (14) for the energy-multiplier ratio together with H — 0 solved 

iteratively for the ratio of X to one of the energy multipliers. 

X 

The segment of the separatrix which is the locus of generalized corners for 
barrier trajectory pairs is simply related to the corresponding segment of the loft- 
ceiling-match curve discussed previously for the fast-evader model. The relating 
property is that tandem -loft trajectory pairs originating on the separatrix segment 


10 


must remain, within a region in energy space characterized by the existence of at 
least one dual - hod ograph- figure intersection and reach the matched-loft-eeiling 
curve, otherwise a transition to 3-D motion and successful evasion is possible. 
The argument is similar to that for the sufficient condition of Ref, 1. In the 
example of aircraft R chasing A , the entire separatrix maps into a single point 
on the loft-ceiling-match curve which Is within the 'region of pursuer-superior 
sustainable turn rate (shaded In Fig. 3) and hence well within the intersecting- 
dual-hod ograph region. 

It appears possible that the locus of generalized corners of barrier trajec- 
tories in energy space may sometimes be a composite including segments of the 
separatrix, the boundary of the intersecting -dual -hodograph region between the 
separatrix and the loft-ceiling-match curve, and the loft-ee Lling-match curve 
itself. Along the intersecting-dual-hodograph boundary, the multiplier ratio 
analogous to (13) and (14) may be obtained from the normal to the boundary, 
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THE COMPUTATION OF BARRIER SURFACES 


Families of trajectory pairs integrated backwards in time numerically define 
the barrier surfaces. All start at Ax~° and at pairs of specific energies along the 
types of curves in the energy space previously discussed. Multiplier values are deter- 
mined as discussed. In the preceding sections. The present section discusses the pro- 
cedure and sequence of calculations employed with computer programs such as those 
of Refs. 4 and 5. 


Perform first a comparison of energy rates and sustainable turn rates with 
energies matched over the energy range for equal loft-ceilings. The computation of 
the barrier surface for the sustalnable-turn-rate-superior aircraft as evader proceeds 
entirely from the boundary of the capture region. If this lies completely above the 
loft-ceiling-match curve, the boundary is defined by matching of maximum instan- 
taneous -turn -rates; if not, some portion of the loft-ceiling -match curve itself may 
form a portion of the boundary. The energy multiplier end values are zero for 
maximum-instantaneous -turn-rate match and given by (13) for loft-ceiling match. 

The value of X^ is found in the latter case from H = 0 . The program of Ref, 4 

does this by a. scan followed by an iteration, if two values of X emerge from this 

X 

computation, one of them may be of interest when the aircraft pursuer/ evader roles 
are reversed. - 


Forthe Sustainable -turn-rate -superior craft as pursuer , segments of the 
loft -ceiling-match curve, the separatrix, and the sustainable-turn-rate boundary 
may form part of the locus of generalized corners of barrier trajectories. If there 
is a region of tandem/loft motion that leads to loft-ceiling match, the separatrix 
should be found, then a composite locus of generalized corners for barrier trajec- 
tories by examining the tandem-loft trajectory family relative to the sustainable- 
turn-rate -superiority region of the pursuer. 


COMPUTATIONAL RESULTS FOR A VERSUS B 

Some results are presented In Figs. 1, 2, and 4 of barrier surfaces for 
turning duels between aircraft A, a hypothetical Mach 3 design, and aircraft B, 
a version of the F-4 used for illustrative computations in Refs. 1 and 2. 

♦ 

Figure 1 applies for A chasing B ; it is the same for both fast-evader and 
neutral-evader models. The upper boundary of the capture set is the instantaneous - 
turn- rate -match curve, while the lower is the loft-ceiling-match curve. Integrations 
of trajectory pairs never develop appreciable positive Ax values except at high 
energies; even in the region most favorable to A, only small angle differences 
can be closed successfully. 

Figure 2 is for fast -evader modelling. B chasing A. The trajectory pairs 
defining the barrier all pass through the lower portion of the loft-ceiling-match 
curve. 


Figure 3 shows the region of tandem motion for B chasing A and the separ- 
atrix arising from bifurcation of a tandem/loft trajectory that divides the region. 

All of the tandem trajectories considered are of tandem/loft type, since all chattering 
singular arc segments are screened out by failure to meet one requirement or another 
(Ref. 2). The pursuer-sustainable-turn-rate region is shown, shaded. 

Figure 4 shows the barrier surface for B chasing A , neutral- evader model- 
ling. The trajectory pairs defining the barrier pass through the separatrix, only a 
small upper portion of the separatrix proving barren. The difference between the. 
barrier results for fast-evader and neutral-evader models is large in this example 
and implies that the tandem /loft pursuit tactic is important. It should be noted that 
this may be untypical because of an anomaly in the otherwise realistic data for 



aircraft A : was taken as unity Independent of Mach number; this is an es- 

raax 

pecially optimistic figure supersonically. This leads to an energy-rate advantage 
for B when both craft fly at m atched-loffc-ce iling altitude, lienee favors the. lofting 
tactic for B as pursuer. 



RESULTS FOR B VERSUS C WITH PARAMETERIZED TTIUUST 


Shown next are some results for aircraft B (F-4) against C , a version of the 
F-5. B can reach substantially higher energies than C (11 QIC vs. 70 K); however, 
it is inferior to C in both maximum sustainable turn rate and maximum instantaneous 
turn rate at matched loft- ceilings. The mismatch in capabilities between B and C 
resembles somewhat the mismatch seen in the preceding example between A and B, 
but is not as extreme. The trajectory families are simpler in two respects: the 
opportunities for successful tandem -lofting pursuit are negligible (and have been 
omitted entirely in the computations), and the relationship between maximum instan- 
taneous turn rates is such that htgh-turn-rate spirals to low energy occur only for 
B chasing C . Le, , C can capture at liigli energy. Taken together, these charac- 
teristics result in simple pairing of trajectories in all cases; there are no multiple- 
subarc trajectories in the families. 


Figure 5 presents barrier-surface results for B pursuer and C evader. 
Successful pursuit occurs only for C in the low energy range and for B with a 
large initial energy advantage in combination with, an angular position advantage. 

Figure 6 shows barrier data for C chasing B . Generally, an energy ad- 
vantage is needed by C to capture B before the chase reaches B f s high-enei-gy 
haven. Given such an advantage, C can close fairly large angular gaps and capture 
when the Initial energies are low; however, at high initial energies, only modest 
closures can be effected even with large energy 7 advantage. 

A computational study of barrier behavior as Bis thrust level is increased 
was carried out. An important guide to the behavior is provided by the sustainable - 
turn-rate boundaries: in energy space depicted in Figure 7, Configuration B2 has 20% 
thrust increase over the basic value over the entire Mach number/altitude range, 


is 



B3 30% » B4 40%. The situation is similar to that of the B/A example as in 
Fig. 3, although simpler in respect to the absence of tandem/loft motion. As the 
thrust is increased from 30% to 40% over basic, B develops an overall 
superiority in maximum sustainable turn rate, energy choice open, and the barrier 
surface for B pursuing C disappears, implying eventual capture in a long- 
duration chase irrespective of initial conditions. Barrier surfaces in Figs. 8, 9, 
and 10 show the shifts in capture capability as B's thrust is increased. The 
results shown are only rough approximations on account of coarseness of mesh 
of the trajectory data. 
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RESULTS FORA TRAPE Z OIDAL- WING . RPV CONFIGURATION' 


Some calculations were carried out Cor the high-performance RPV configura- 
tion of Ref. 6. Maximum instantaneous turn rate, maximum sustainable turn rate, 
and maximum loft-celling rate versus energy are shown for this configuration, 
designated D, in Fig. 11, where the corresponding characteristics of aircraft B 
are also shown for comparison. The superiority of D over B is overwhelming; 
the sufficient condition of Ref. 1 is met and eventual capture in role-determined 
D/B engagements is insured irrespective of initial conditions. 

A similar comparison of liodograph data for D and A is shown in Fig. 12. 

D is seen to be superior in the dogfighting energy range; the superiority wanes 
above 100 K energy because of inlet off-design characteristics. The barrier re- 
sults of Fig. 13 show A incapable of capture except at extremely low evader energy 
and then only with tremendous energy and angular position advantages. The barrier 
surface of Fig. 14 indicates high effectiveness of D pursuing A except at high 
energies near A's haven. Fragmentary trajectory data in the energy range 
60-100 K indicate that the poorly -defined barrier surface may actually fold over 
into the region below the loft-ceiling-match curve {marked No Capture ) for large 
angular separation; there is some difficulty in representing the surface graphically 
here, in addition to a coarseness-of-mesli problem. 
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CONCLUDING REMARKS 


A procedure Las been presented for the calculation of barrier surfaces for 
differential-turning duels and examples given. Three computer programs, each 
fairly complex, were used and a certain amount of crossplotting of intermediate 
data was required in the first computational, example; however, considerable 
streamlining of the programs and computations was carried out in the process of 
the subsequent applications and the generation of barrier surface data can be re- 
garded as a. practical proposition for engineering applications work. 
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ABSTRACT 

This manual is a guide for using the Aircraft Energy-Turn and Tandem 
Motion programs, digital computer programs for use in differential-turn com- 
putations. The program is written in FORTRAN IV for use on a CDC 6600 
computer. With slight modification, the program may be run on any digital 
computer which will accept FORTRAN IV. 

The program inputs and outputs are described herein and examples pro- 
vided. in addition, the program options are discussed and overall flow charts 
are presented. 

The programs described were developed partially under U. S. government 
contracts: F44620-71-C-0123 with the Asst. Chief of Staff, Studies and Analysis, 
Headquarters USAF; N00014-73-C-0328 with the Vehicle Warfare Office, Office 
of Naval Research; and NAS 2-8738 with NASA Ames Research Center. 
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INTRODUCTION 


The Aircraft Energy Turn program is an outgrowth of a digital computer 
program which was designed to go parameter studies of variable -altitude turns 
obtained by numerical integration in reduced-order ("energy’*) approximation for 
supersonic aircraft, including the effects of constraints on altitude, dynamic 
pressure, Mach number, lift coefficient, and normal load factor. The general 
theory of reduced-order approximation for aircraft maneuver optimization is 
described in Ret 1. The program described herein provides for the simultaneous 
Integration of two aircraft trajectories, uncoupled except for control of print-out 
and run termination. The present user’s manual overlaps that of Ref. 6, which 
describes a similar single-aircraft program. It dt scribes, in addition, a re- 
lated "Tandem-Motion" computer program. The two programs are used in con- 
junction with a third, "Hodograph", computer program, described in Ref. 7, for 
differential-turn computations. 


The following tables are utilized in the programs; vs. Mach, 


vs. Mach, rho vs. altitude, vs. Mach, speed of sound vs. altitude. Thrust 


tables, which are functions of Mach number and altitude, are also used. Internally, 
the tables are fit by cubic splines or cubic spline lattices (Ref. 5), with all "mo- 
ments" (second derivatives at mesh points) stored In tabular form to be used for 
interpolations during trajectory calculations. The desirable feature of spllning 
for the present application is enhanced smoothness (continuity of second derivatives 
is maintained), which facilitates the generation of families of solutions of the state/ 
Euler differential equations. The speed brake drag coefficient, C L » is introduced 
as the product of a factor and the zero-lift drag coefficient. A preprocessor Is built 
into the programs so that the user may change any of the tables with input. This is 
explained in detail in the Input section. 
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TECHNICAL SUMMARY. ENERGY-TURN PROGRAM 

This section summarizes some of the main features of the analytical material 
presented in Ref. 1. The state variables in the basic system of differential equa- 


tions are the specific energy, E { = h+v /2g), and heading angle, x* The roles of 
control variables are taken on by altitude, h , bank angle, p, and throttle setting, 
O^r^l. The Euler system Is given by 



The Hamiltonian, H , is given by 


HE x E XM + x x f tan , 

The augmented Hamiltonian, H , is defined: 

5 

>**.. - •. S = ? + Z\^i 

' 1=1 • 

where $ are the Inequality constraints: 


1! 

h-h T ^o 

terrain limit 

"CD 

ro 

II 

q - gp (E-h) a 0 

dynamic pressure limit 

*3 “ 

M - M* 0 

Mach number limit 

II 

cos , 1-^*0 

C, limit 

'‘max 

K * 

COS fj, - cos ft a 0 

normal load factor limit 
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Further details concerning the determination of the multipliers X^ and the control 
variables are given in Ret 1. 


It is instructive to consider two basic types of maneuver which are limiting 
cases. These eases are: (1) symmetric flight. In which the final energy of the air- 
craft is to be maximized, and (2) hard turning flight, in which the final heading angle 
is to be maximized regardless of energy loss. The Hamiltonian is a weighted sum of 

energy and heading rates, namely: X E and X x • The program finds the altitude, 

S. X 

h , at which the Hamiltonian is a minimum by searching versus h at constant energy 
level. 


The first case, symmetric flight, can be generated by making X negative 

A 


(say, -1. 0) and X zero or small negative. Because x Is generally much 
.X max 

smaller than E for a given value of E , X ~ -1. 0 is considered a small value, 
max X 

(Actually the multipliers X and X need be scaled such that H = -1.0, but since 

“ X 

only the ratio affects the computations, this scaling can be postponed.) A typical tra- 
jectory might start with a low energy level, say E = 5 to 15,000 feet. With such 

values, the aircraft executes a generally accelerating, climbing path beginning at 

... #' 

sea level, choosing an altitude for which E is maximized. Two typical symmetric- 

flight trajectories are shown in Fig. 1 in an altitude/veloeity chart. The trajectories 

are shown superimposed on the velocity-altitude envelopes corresponding to constant 

values of E . Constant energy lines are also shown in the background as in the pro- 

* # 

cedures of Refs. 2 and 3. The aircraft climbs along a path of maximum E , (bE/bh)„~0 . 


There is a subsonic and a supersonic branch of the function (bE/bh) =0 , which overlap 

* - - , JJ# 

in the energy range of approximately 32, 500* < E< 45, 000 t ; i. e. , there are two minima. 


The branch taken depends upon how the search for minimum of H is initiated; in other 
words, the program does not have the capability of locating a global minimum. The 
times at which the points along the trajectory are reached, as well as all values of the 
variables and airplane parameters, are printed at each time-step by the program. See 
the section labelled Output for further details. 
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20 Square* to the Inch 
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The program does not calculate the altitude-velocity envelope* energy -rate 
and turn-rate peaks; this is available from the companion "hodograph” program of 
Ref. 7. Figure 2 presents the peak values of energy rate versus E for the sub- 
ject aircraft and illustrates the subsonic (first) peak and the supersonic (second) 
peak. The overlap range is, as noted previously, between E- 32,500 and 45,000. 
The first trajectory (lower) of Figure 1 transferred from the subsonic to the super- 
sonic branch of constant energy level E ^ 34, 000 , while the second transferred at 
E ss 45, 000 . The difference in these eases depended upon how the initial altitude 

iterate was introduced into the I -D search; selected by ITS=0 , where h is 

guess 

the minimum of h, , h n , and h_ , or by ITS-1 , where h is the value h cal- 
12’ 3 guess 

culated in the previous time interval. The second trajectory (upper) required a longer 
time interval to reach the final energy state than the first, since the first peak was 
followed for too long. Ideally, the transfer point would occur at the intersection of 
the first and second peak curves, E a? 36, 000 feet. 


The second case of interest is that In which maximum turning rate is sought 
irrespective of energy rate. In these cases, the Hamiltonian consists of the term 


X £ tan ji. Trajectories having this character may be generated by taking X as a 
X » n X 

very large negative number, of the magnitude: X = -1x10 , 

X 


where n & 8. The Hamil- 


tonian is minimized when the ratio (tan ji)/V is a maximum, and this corresponds to 
the so-called corner velocity described in Ref. 1, where the load factor as limited by 


Is equal to the maximum allowable load factor of the airframe structure. The 


max 


corner-velocity locus is not normally generated by the program during a trajectory 
calculation, but an option is available for executing this calculation. 


Figure 1 shows a "corner-velocity” locus in turning flight in the Mach number- 
altitude plane. As shown in the figure, in the region to the left of the corner velocity, 
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20 Square* rq the Inch 
















the maneuvering capability of the airplane is limited by , while to the rights 

the structural load factor limits maneuvering. maX 


The general optimal turning maneuver is Intermediate between the symmetric 
energy climb ease and the ,, corner-veloctty ,t case; however, there is more complexity 
in the general combined turning case than in either of the limiting cases. Thus the 
Hamiltonian function may exhibit three relative minima versus altitude h at constant 
energy, and the transition behavior may be quite intricate. Because h, in its ap- 
pearance as a control variable, is not "the real thing , n having arrived by way of an 
order-reduction approximation, the analyst using the computer program to generate 
families of Euler solutions may be interested not only in the global minimum of H{b) 
but in the relative minima as well. Trajectories so generated correspond to so-called 
"strong" and "weak" extrema, respectively, and complicate the survey of families. 
Transition between minima should be carefully monitored and controlled by the user 
constructing families. 

The user may elect a mode of operation that provides the global minimum of 
H(h) and hence "strong" extrema. 1 2 3 This is accomplished by a preliminary scan of 
H values at 500 ft altitude increments over the complete altitude range. While it is 
computationally expensive, the alternative of controlling transitions off-line is laborious 
and slow. 

Where "weak" minima are of interest, a measure of control over which mini- 
mum of the Hamiltonian the program selects is provided the user by two optional input 

choices described in the next section. One of these permits choice of first guess for 

2 

the one-dimensional search as either the value of h resulting from the previous 

1-D search at the preceding integration step or a value slightly below the minimum 

3 

altitude limited by terrain, dynamic pressure, or Mach number. A second option 


1 ITS = 3 

2 ITS = 1 

3 ITS - 0 
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consists of a choice between normal search logic described In Ref. 1, on the one hand, 

4 

and a logic choosing the minimum altitude whenever the switching function X> 0 

d ; H ** 

(cut throttle) and the slope ^ given by Eq. (123) of Ref. 1 is positive, this restricted 

(and possibly "weak" minimum) choice reverting to normal logic should either sign 

change. . ; 

A mechanized search for the multiplier value X yielding a specified energy 

X 

at a fixed final time Is available upon input option. The iterative scheme employed 
is minimization of the square of the difference in actual and specified final energies 
via one -dimensional search using golden section and cubic fit techniques of Ref. 4. 

The program provides for simultaneous integration of energy and heading- 
angle histories of two aircraft for convenience in differential-turn computations 
(Refs. 7 and 8). The computations are independent; however, the prtnt-out interval 
is controlled in part by turn-angle difference and the integration is terminated on 
multiple criteria which Include violation of prescribed bounds on the specific energy 
of each craft. The choice between integrating trajectories for one or two aircraft 
is controlled by input, as discussed in the next section, as is print-out interval 
. and run termination. — - - 

For the determination of X initial values that make the Hamiltonian for 

X , 

the two-vehicle system vanish, hence result in grazing (barrier) trajectories, an 

- iteration procedure is supplied. This minimizes the square of the Hamiltonian • 

via golden-section search. 




NAMELIST/MINP/ATAB, HATB, CDOT, XMOT, CDCLT, XMCLT, 
RHT, HRT, TCLH, YMCL, TRXH, XMTRT, HTRT, NMA, NCDO, 
NCDCL, NRH, NCB, NXT, NHT, IDEN, EPSCF, EPSGS, LT 


EDEN is dimensioned 6 and each slot is associated with a table ID. The 
tables are as follows; 

Number 

ED Table Dimension of Points Description 


1 

ATAB 

20 

1 

HATB 

20 

2 

CDOT 

20 

2 

XMOT 

20 

3 

CDCLT 

20 

3 

XMCLT 

20 

4 

RHT 

20 

4 

HRT 

20 

5 

TCLH 

20 

5 

- YMCL 

20 

6 

TRXH 

20,20 

6 

XMTRT 

20 

6 

HTRT 

20 


NMA 

NMA 

NCDO 

NCDO 

NCDCL 

NCDCL 

NRH 

NRH 

NCB 

NCB 

NXT, NHT 

NXT 

NHT 


a (speed of sound 
h (altitude) 


S 

Mach No, 



Mach no. 
rho (density) 
h 


Maeh no. 


TR (thrust) | 


Mach nos. = 
Altitudes = 


Maeh no. 



cols. 

rows 


If IDEN(7) “ 0 , the CL table is read in directly. 

C2 

If IDEN{7) ^0, the table is read into the CDCLT table and the 
program inverts it before preprocessing. 

If all the tables are to be preprocessed, IDEN(I) £ 0 , 1 = 1,6. 

If the user expects to use the same table or tables for many runs, he may 
want to avoid the preprocessing by reading in IDEN(I) = 2 on the first run. 
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This will give a printout of the moment table or tables, which may be put 
Into data blocks in the succeeding runs and then EDEN (I) should be set to 0 . 




LT = 0 Indicates end of run 
LT = 1 First aircraft tables are preprocessed 
LT - 2 Second aircraft tables are preprocessed 

LT =■ 3 Iterate for minimum HI< X ), Ht{ X ) = H 9 ( X ) - H ( X ) - H* 

See PROG1; Summary of Subroutines 
LT = 4 Integrate two-aircraft trajectory 

LT = 5 Integrate one-vehicle trajectory only _ 

EPSCF Epsilon In ONEDIM to measure closeness to origin 
EPSGS . Epsilon for convergence tests in above . . . 


B. Main Inputs 
1. Aircraft One 

NAM ELIST/MINPU T /GO, H, S, W, E, X, XLE, XLX, TMAX, DELT t EMIN, 
EB3, BDLT, EPSGS, EPSCF, NEQ, NPRN, IND, INTB, QBAR, QMAX, ICB, 


DLLME, 

HSTR, NC, XBAR, XMBAR, XNB, CNCLC, ITS, INCR, HT 



Notation Explanation 

Units 

GO 

Gravitational constant (32. 16) 

ft/see 2 

H 

Guess for Initial altitude . / . .U 

ft 

S 

Wing planform reference area 

sq ft 

W 

Weight of aircraft 

lbs 

E 

Initial energy , 

ft ; 

X 

Initial x angle 

rad 

XLE 

Initial X„ 
E 


XLX 

Initial X 



X 


TMAX 

Maximum time of flight 

sec 

DELT 

Integration interval 

sec 

EMIN 

Trajectory terminates If E a EMIN 

ft 


ORIGINAL PAGE IS 
OF POOR QUAIOT 


EB3 

HDLT 

EPSGS 

EPSCF 

NEQ 

NPRN 

JND 

INTB 


QBAR 

QMAX 

ICB 

DLLME 


HSTR 


If E < EB3 , bypass A calculation ft 

u 

Initial altitude step size for the regular ONEDIM routine ft 

Epsilon for convergence in the golden section of the 
ONEDIM routine 

Epsilon for closeness to origin 

Number of equations to be integrated 


Number of integration steps per printout 

Print option Indicator (to be explained in the Output 
section) 

0 - last case was processed 


1 = always calculated 

2 = B Ae not calculated when X > 0 and bH/fch > 0 

45 Jci 

Q boundary for calculation 


Q boundary for cubic fit in ONEDIM routine 
If Q > QMAX , use cubic fit 
If Q £ QMAX , use golden section only 


0 - 




1 - Gj^~ TCX.H (Mach) 


First step in ONEDIM on ®( X ) 

X 

If initial X = 0 , program will start with -DLLME and 

X 

terminate if it cannot find a minimum 


If initial X < 0 , program will reverse directions if It 

X 

cannot find a minimum 
= H* 


NC 

XMBAR 

XNB 

CNCLC 

ITS 


Number of samples allowed in INEDIM on HI( X ) 

X 

Mach boundary for ^ calculation 


n (maximum structural load factor) 
Speed brake factor, where C^= CNCLC 


% 


First guess for altitude in 1-D search 

0 - taken from previous ^ , or ^ solution 

1 - taken from previous 1-D search 


500. 


2*E 


3 - scans an altitude range for minimum Hamiltonian 

ORIGINAL PAGE! IS. 
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INCH 


HT 

XBAR 


Number of steps in scan mode 
Maximum altitude allowed 



Xg ^ XBAR , trajectory terminates 


ft. 

rad 


2. Aircraft Two . 

NAMELIST/MINP2/H, X, W t E, XLE, XLX, EB3, HDLT, XNB, INTB, 
QBAR, QMAX, ICB, XMBAR, CNCLC, ITS, INCH, HT, S, IND, EMIN 


See above for explanation. 


C. ERROR CODES 
Printout 

Number Location 

40 Table look-up 

100 Scan mode: lower limit exceeds upper limit 

111 \ Aircraft 1: h>E 

112 Aircraft 1: cop p* > 1 

121 Aircraft 2: h >E 

122 - Aircraft 2: cos ji* > 1 

334 • Subroutine SPINT: Interpolation for X ~ 0 
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Sample Input 


Two-Vehicle trajectory: Aircraft One 

Aircraft Two 


SMINP 

NMA*20, 


( i EHi O t t 1 * t»i iU-.t i 

ncoo-19, 

CiJOr = . O20%,,O2C)5 f ,»2Q5, . Q.?H 5 , , C S * . 0213 , . 02 11 , . 0242, . 0 J24 , . 0 15 9, . J 374 , . 0394 . 
.i>3d8,«'336p,,Q3dy,*Q397f*j9B3,»39Jif.ii.9 < 3 3f0», 

X H0T= 0, , • 2 , « 65 , * 75 , • 60 , • tip , , 9 0 , • QG , 1 » 00 , 1 • B5 , ! . 10,1.20,1.40, 1.60,1. SO, 2., 

2<2|2t4i?«6iQ»f 

NC-B= 14, 

XMCLT— 0 « t «2i *85, ,75, *30, • *% , .9 3, ,95,1*0, i . 05, }. . 10,1,20,1, 90, 1,60, l, 10 , 2 • , 

2 . 2,2*9, 2*8*3. , 

NCDGL-19, 

GuCl T~ * 198, * 19 ft , , 194, ,2 ill, . 2 (3 5 , , 211, ,210, ■ 22 3 , * 23 j , * 2% 2, *265,, 297, . J63, . 4 5 S , 

• 544 , *-64 , . 7u 4 , .7 44, , 744 , 3 . , 

Y HCL = 0 , , • 1 , ^3 , . 5 , • 7 , * 3 , 1 • t , l . 5 , i * % , 1 , 7 , 1 , -j , 2 . 1 , 2 • 3, 2 • 8 , 3 * , D . , 0 , , 0 • , 0 * , 3 . , 

TCL4-.38, *3-3, .31, • -9 1 , • 7 4 , *93,. 90 , * 75 , • 54 ■< * 45 , . 4-J , • 3 5 »-* 3 B , » -3, 

N X T = 14 , ;jH f - 1 0 , • ' ' ... 

X9T*<T = 0» ,, 2, ♦ 4, 3,1,, 1.2, 1.4, 1.8,1. 3,2*, 2.7, ■* * 4 , ? * 8 , 0 • , 3 • , 

HTRT^tIOOOOO. , Q * flOtiOG* , 2 y 3 0 0* , 3 2 y U 0 . , 4 tj 0 U 0 * , 5 0 <J U 0 * ,60000. ,7 0003. , 

200 0 JO . » 0. • y • , 0, , 3 . , D . , 3 . , 0 . , 0 , , 0 . , 0, , 

T " X H 1 1 , l ) - 1 2 0 2 a 0 . , 1 2 0 2 0 0 .,1174 0 0 . , ID A6Q3 , , 119? 3 5 . , l !i 20 00 , , 33 4 ll 0 • , t ’ C 0 .. V 6d4 . . 


; u d U . , i. i t*,u. , U . * j t , . , y , 

7 X .9 ( 1 , 3 ) = •? i B 0 . , 2 8 3 J • *29 0 - • , 3 ] GO . ,383 > . , 49 00, ,560*?, • 70 OB , , 3 4 0 U • , 33-10, , 
u&OQ*,lU'du*j*,1320 3*,78'3Q«,‘j»,v»,!i«,G.,J»,Q., 

^XH 1 1 , '3 ) - 2 0 00 • , 20 33. , 2 3 0!' . , 27 Q 0 * , 24UU . , 2 **. QfJ , , j >o o 4 tiOG • , 4S 30 , ,550 3 , , 
200 . ,5403 . , fi 0 OQ • ,45 0 J . ,J.,0.,G. , • , C « , 0 . , 
r-XN (1,13) =?G*0. , 


V*** — 

OF POOR 
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SMIMP 
L T = 1 , 

NMA-2Q , . .... - ^ * 

At A0-i49'fi . *11 16 • Via 97 * , 1Q77V, 1C57. *11337. *1 016 . ,994 .7,963.1 ,96-3. 1,96 8.1, 

913. 1,968. 1,968.1, 96ft. 1, 958. i ,968.1,984.2,1004. ,1C56. 43, 

HAT3 '-1 » E5 , 9 « ,5 . E3 t 1 *.E4 , 1. 5£4* 2 » £4, 2. 52 4*3. £4, 3.6432 6, 4. E4 ,4.5E4*5.E4, 
5.SE4,6*E4,7.E4,8.E4,a.Ea2Eiir f 9.E4- r i-.E-5,2.FS f - 
N’H=2u, 

RHT = 6321 • E- 5 ,2 377 •£— 6 ,20 48 .E— S ,1755 .£-& , 149S.E— 5,i266.£" , &,10&5. E“ 6 , 

8893. E-7-, 7Q61. E-7-, 5 851. &-7,46Q 1 . 6- 7 ,361 ft , E-7 , 2 345 . E- 7 , 22 38 .£>-7 , 

1334. E-7,8556.E-a, 7764. £-8,5151 .£-8, 3131. t-8, 94 -t-ft, 

BRT— •* 1 ♦ £5 , 0 . , 17.13,1 • E 4, 1 . 5E4 , 2 . E4, 2.5E4, 3 * E4, 3 • 60 9E4 ,4 * E4-, 4. 5E4 , 5* E4 ,* 

5 . 5E 4, 6 . E4 , 7 . £4 , 8 VE4 , 3 . 2(|2E 4 , 9 . E4 , 1 . £5 , 2 . E5 , ' ‘ 

hC3= 16 » 

Yf’!GL-0*,.2,.4,«5,.7,*3,.9*.9 t j,l.,1.05,l.l,l*2,l*4,i.6,1.3|3.,4 + Q., 

T CL H = *86, ..36 , .8 6 ,..79, . 76, .74, . ft 4 ,1 • Q4-, 1 » 2 » 1 • 2 7 , 1. 25 ,1.13, .85*. 65 , #-5 , . 5 , 

4*0.,. . .. ..... 

NC90=16, 

X.mOT-O., .2, .4, .6, .7 , .3, . 9, • 9? * 1* ,1.05,1* 1,1* 2,1 • 4,1 ■ 6 ,1 *3 , 3 • * <4* Q • - 

CUQT- • 02, .02,. 0195, .0192, .0183, . 0189, . 02 03, . 0 27,. 0 38, . ff#OS 0 406, . 0409 , 
.0414,, 0413, .04U,. 0413,4*0. , .. 

ncocl=16, „ 

X HC L T — G • , .2, * 4 , . 6 , .7 , .ft, .9, .9^ ,1 . ,1.06,1. 1, 1.2, 1*4, 1.6, 1.3, 3. ,4*0., 

C JCL T— 5 *4, 6.4, 6.4, 6*4, 6.6, 7*6, 7. 1,6. 6, 5. 5, 6.4, 5*3, 4*6, 3* ft, 3. ,2*5, 2*5, 4*0 . , 

UXT = lfl,MHT = 9,- . 

Xi'i T RT .,*2,*4,, 6, .8,1. ,1*2, 1.4, 1,6, 1.8,13*0., ~:* 

BUT =^15 30 0 , , 0 ,, 15 0 0 0. , 25000 .,360 ft 9. ,40000. ,450 00 . ,5 0 U 00 . , 1750 00 . , 11*0 . , 

T !i X-H 9 1 &O-.^9i.0-Q . ,-9901. , 1 0 70 0 - Vi 1 50 0., 1210 0. * 132 00. ,14300. ,15430 . V- • 

1653 0 • , 10 * 0 . * ' 

r 2 X H t i , 2) = o60Q,, 6600 * , 72 03., 73 uJ.,87Qu», 9500., 1030 0«, 11100., 11900., 

■ -12-70 ;; . r£s r~. -:\ ■■ — r ,_ . . 

TRXH (I ,3)=4100. , 4100 . ,4 500. , 51 0 0 . ,5900 . , 63 00 . ,740 0 . , 790G . , 84 00 . , 390 0 . , 

10*0., - .. ..... ... 

TRXli 1 ,4 )— 2-50 0 .-,2 501 *^,-3-0-00 . , 35 0 0. ,4100* ,50 00 ,,o 3 0.0 * ,6900. ,7 ft 00- .-,87 QO * -- - -~ ; 
10*0., 

TRXH (1,5)-1UUU., 10 00. ',16 on. ,2200 .,2600., 3300 . , 4 10 C., 4800., 5400. ,5930., 

10*0., - ..... ... 

TRXH ( 1, ?) ~7G 0 . ,700.,, 1200. ,1700 * , 22 0 0., 2303*, 340 *1. *4000., 4600. ,5200. ,10*0., 
TRXH( 1,7) -8 00. ,800. ,1100. ,1400 . , 17QG, , 21 03 . , 270 0 . ,32 0 U . , 37 00 . , 4200 . ,10*3. , ■ 
TRXH (1, ^>=600-., 600-., 700 .-,1000 . , 1 300 160 fl 200 3. , 24GTT. , 2 80 Q. ,3200. ,10*0. , ; 

TRXH (1,91—20 *G . , 

I Uc S=l, 1,1, 1,1, 1,1, 

IE HD •• 
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$; M iI N P 

tPSF=i . g-5 , t-P5G3 = l . E-5 1 
LT = 4, 

$CU3 
SHIN PUT 

GO- 32* 2»H =13330. ,3=170. , U—.JflQU* , ' =40003 ,,XL£=^i. ,XLX=-t*, 
TMAX=203.»DCLT=-1. , 

CeAR=4QQ9 . , IC5=Q ,XN3=5. 7 391 ,X 10A:i=-3, ,CNCLC=1. 5, r :*T3 = 1 ,E13 = 15 lj]Q. ( 
ITS=3,HJLT = 10 . ,QMAX=17Q 3QO . , t NCR=5 00 , EMI M~Q 0 Q 0 1 . R , Ei’SGS- • IE- i % 
W=55 4 £ . , XN 3= 7 ♦ 3 3 , Q3 A R- i * 0 ■:! , , X 1 11 A >1=1.72, 

£PSCF=.iE-l *, , IN'JsltNPnN = 5,. ! -JEQ“b ,h \~ 2 , , 

XbAR=S., 

XLX=-1. , 

E=lfJ5J3 .,XLF=-.1E-7,XLX = -U , 

$EMO - 

SMI 'if *2 

S=5 3 0* , = 37 52 3 • » X = fl . , IN0R=5 00 , CMCLC=1. ,U'3A 3 =2l5 7, V y , x 4 4 AR=2 . 

XN-B- 7 • ,XC3=1*H=25U00. , 1 1 3=3 , 1 N T 0=1 , HOLT = 5 0 • , 

£03 = 130 000. , 

XLE=-2. ?1G5Q 7 7501 43147, E = &7395, 7 55 . 4S6957,4MAX=2m\ , 

XLX = '. 47555 32 015 72544 5, = 120 0 0:. , 

£=75 0 0.0,XL£=».lE-7,XLX=-t . , 

I MQ= 1 , 

2»EHG 

sminput 

I Hf 3=3 , 

SENG 
SMI ii P 
LT = 0 , 

3 EMU 


SSS'wasw- 
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Sample Input for a two-vehlele Iteration followed by a trajectory: 


ZHV*P 
Li = 2 , 

ATA0=149o. , ill 


1697., 1277., t 
9 fc b • t * 3 S 'i . 1 » 9 o 8. 1 * 9 & 0 * 1 *9SS» 1 * 
Hfif3 = -l.E5t‘^ • * 5 . c 3 , 1 . E4 * 1 *5 £ 4 , 
5 • 5 i 4 i 3 • E 4 9 7 » E 4 1 3 • £.4 , 3 • 2 p 2E4 * 9 

NX 8 = 20 * . 4 .;.: • • a... . .. 

RHT=SS2i.E-3*2377.E4&,2043.E-5 
*398. £-7, 73 si. £-7, 5351. 2-7,460 
1334»E-7*<i5 ! 36.2-3, 7 764*5— -it 5 1~ - 
H^T=-1.-E6, 0.,5.£3,1.E4 


0 -> 7 * *1037*5 1Q1 & » i 99 4 *7 *933*15 9o 8*1*9 
93b, i*9 8 4*2*1004**1 Opb* 4 3* 

Z . £4 ,?* 5£4,3. E4, 3.609E4, 4. E4 ,4.524,5 
• ii4*l*E5,2.25, 


9 8 . 1 , 
• £4 , 


Pa I * l iiO* ji * ; 1 * ?E H* - 

5.5*4, 5. E4, 7. £4, 3. 64^.» .iG2E4* 3 
NC03-19,- • - 

C 0 0 T = . 0 2 S 5 * . G 2 95 , .02 05* .02G5, * 
.C335* *j385, . ■] 337 » • 0 397 5 • *040 i, 
-X.H-OX — 0 * 5 * 7 * * ,5 7s)-*- * - 8 G 5 *-35-*-*-9 > 

2»2,2*9,2.q*U»* 

NC3= 14* - ■ _ '• 

X hCl T = 0 . * * 2 » • 6> >* . 75 , * 3G , * •■‘•5 * • 9 
2 , 2 * 2 . t *, 2 » 6 » 0 *» . ■ 
mco:l = i3, ■ 

C9GI T = . 198 , . 1 m 5 . li^5 *.2G 1 5^2 0 5 
•54 4 *.54**744, *744*. *744 * 0 » » 

V 8- • 5’?0.*.1,*3,.?*«7**9*1.1*1. > 
T 'j-L > i — * S 3 * » 3 :► *—* -8 3 , 31**735 *5 8 * • 


*l785*g-otl496.g-a,1266*e-6,1065.E-5 
1 . £>7,3613. F -7, 2 345. £-7 * 22 33. £- 7 + 
i. £-8, 3131.2^3, 94. E-3, - - 

£452. £4 * 2*5e 4 « 3, l£4» 3 • 50 924 * 4. £4 * 4. 5E4*5* 
. E4, 1.E5* 2.E5, • ■ - • 


0 2^5** 0 213 ,.0 210, .0242*. 0324*. 0 359, . 
* j 4 0 3 * *040 3* 0 . * 

* « 95 * 1 ,00 * 1 . 05 * 1* 10,1, 20* 1*40* 1* 50 * t 


£ 4 , 

3374**0334* 
* SO, 2 * , ■ 


U ,. 45,1,0,1.05*1. 10,1. 20,1.40,1.^0,1.30,2.5 

, . 211* .213, .223*. 239**252, .265 ,-297,. 353, .453* 

*4 *5, 1.7,1. 9 *2. 1,2.3, 2.5*0. *Q.*0.*Q.,[J,*0.» 

90, , 75, » >4, . 45 , . 4 0* .35, ,30 ,.3* 


i* • ♦ U 


. * 3 . * 3 . * 0 . *- * * 


N>;r = 14 *'IHT = i!}, ,. : 

XHT r = 0. ,5-2 ,-54,. o*. 3 .1. *1*2*1, 4* 1:5-5, 1.-3 , 2, * 2 . 2 , 2 . 4, 2 . 5*0. , 0,- , li * , 0 * 
HT^T = -I0a000 ,* 0. ,103 09. , 2900 0 . , 3 s fife ,40 GO 0,4 5U 3 0 0 . , & 0 00 Q . *70 900., 

200000 .,u.,-v., 3 . , 3 , , 0.5 0 *, 0 . * 0 * *G» , 0 . * 

T 2 X H ( 1 , 1 > = 1 2 Q 2 9 0-54:120 2 0 0- -544-74 0 U-. *1038 0 0. , 1192 0 3 -,10 2000^,-334 0 0:*., 2 


, 0 5 , 0 . 


TRXH (1 
4230 0 * 
TRXH ( 1 
44300. 
T2Xrt(l 
37200. 
TEXH (1 
26490. 
T3XH <1 
1783 0 . 
19X8 (1 
iiia 3 0 * 
T?XH (1 
62 0 0 . ,5 
T * X 8 (1 
ID£N= 1 
SEND 


1 0 2-0 00 ^,- 33400 ^, 2 200 V * 5 i 
41400. ,40200. ,39800. 
•42200 . *4403 0 . , 4300 0 . 



5 * 


1 * 1 * I * 1 , 1 * 0 , 
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ys ! ‘ 


SfllMP 
Lr=i » 
hHA=20, 

A T Ail — 1 4 9b • , 1 1 1 6 * , l 0 9 7 , , 1 07 7 • ♦ t 0 1 ( ♦ » 1 Q 3 7 * , 1 G 1 6 * , 10^ * ? t^B8< 1 t96=i* l i 958*1 1 
9&S. 1,96 ') * 1 ,95*. 1,963.1 ,9v« . 1 ,9a 6.1,9 84. 2, 1 104. ,14 56 .43, 

H A T b — — i*E5,9.,:?*EJ,l.E4, 1 » 5 F 4 , 2 • E 4 » 2 « 5;. 4 , 3 » L 4 » 3 . 6 b 9 if 4 » 4 » E 4 , 4 * 5 E 4 , 3 • E 4 » 
?3.5a4,o»E t *-,7«£4,8.E4,8»<?U : ‘~*A,-5»E4,l.if 5 *-2 • f 5 , 

IKH-20, 

KHT = 63 2l,E-5,2 37 7 .E-6,20 4?’ .E-5 ,1 755.6-6, 1495,F- , t l2f,5.E-&, IQ65.E-6 , 

8393 .E^ 7 , 7061 • E— 7 * 5351. E-7, 46 U 1 . E-7.,3618 .£-7 ,2 8 45 .F-7 , 2?38 ,£>7 f 
1334.1-7, 8556. E-3, 7764. E-8, 5151. t-3, 3131 #E-8, 94. E-8, 

^ “ a » '_6 , D . » L* . -f. 3 , 1 , E4 , 1 • -E 4 * 2 • E 4 » 2 * y c.4 , 3 , E 4 » > * 40 v> l. 4 , 4 » k 4 , 4 * 5 C 4 , 5 ■ £ 4 , 

5*5E4,b.£4i 7 . E 4 » 3 » E 4 , 8 . 2 0 2 E 4 * 9 • E 4 . 1 • E 5 , 2 .’75, 

NC3- 16 » 

T k. L H - . 8 4 , t 3- i • 9 a j .79* .7 a, *74, • 84 ,1.0 4,1. 2,1.2/ ,1.28 , I .13, ,35,. 4 5 , * 5 , .5, 
4*0. , 

MC JO = 14, 

Xit»T — Q., .2, *4, * o * .7 . • 8 , .3, .93, i* , L ,0 b ,1 « 1, 1.2,1 .4, 1*6, 1.8, 3. * V* 0 » , 

0 tj 0 T = . O 2 , . Q 2 , . 0 1 9 ■> , • Q 1 '9 2 , . 0 1 8 8 , . 0 1 1 9 , * 0 2 0 9 , « 0 2 7 , . ■ * ,5 3 * • 0 4 D 3 , . Q 4 *1 6 , . p 4 Q 9 , 
.0414, .0413, .0413, .0413, 4*8. , 

HGOJL=16, 

X'lijL T- 0., .2,. 4, .6, • / ,.3, * 9 , * *) 3,1,, 1*3 5, 1.1, 1.2, 1*4, 1*5, 1.6,3.,'* *9., 

6 j CL T = 6 • 4 , 6 * 4 , a , 4 , 6 • 4 , 6 » b , 7 ,6,7 • 1,6. 6, 6.7, 6 .4,3 *3, 4. 6, 3.3, 3. ,2. 5, 2, 5, 4*0* 
NXT = 10 , --4H ?“9 , 

XMT R T= p . ,.2, . 4 , * b » , 3 , 1 . , 1 . 2 , 1 . 4 » 1 . b , i . 8 , 1 U * 0 . , 

HT.^T =-15900 ., u 15GO0 ., 25000 ., 3. bu 39. ,43300 ., 450 00 . ,50000 , ,175000. .11*3, , 
TkXH < 1 , 1 » =9 1 0 0 . , 9 1 00 . , 94 00 . , 1 ; 1 7 3 G . , li 50G . , 1 2 1 0 0 . , 132 00 ., 1430 9. , 15 4 0 Q . , 
laSOQ. ,10*0 . , 

t*;X.H < i , 7 > =66 O J . , 6600. , 7200. , 7 >0 0 . ,1700. ,9503. , l 0303 , ,11100 . ,11 ICO . * 

12703. , i!i*G. , ■ . . . 

TFXH 11, 3) =410 0 • V410U . ,45 00. , 510 0 ,7590 0, , 59 30', , 7400.77900. ,840 0. ,890 0. , 
10*0.-, 

T^XH (1,4) =2500 , , 2503 . ,30 0-0 . , 3 5 0 Q , , 4 1 0 Q . , 5 0 0 0 , , 6 0 0 0 • , 690 0 • >78 0 ti . ,8730 . , 
10*0 . , 

T^XH (1,5 > =1 a 00 . , 10 3Q . , 1500. , 2? GO . ,2600. , 3333 . , 4 100 ,, 46(10 . ,540 0 . ,5900. , 
10*3 . , 

Tr^XH (1,5) =7 30, ,700. ,1200. ,170 0 . ,2 CO 0 , , 23 G 0 . , 3 4 0 j . , 4 0 0 0. ,46 00 . ,5200. ,10*0, 
T3XH(1,7)=3U0. ,333.41103. ,1400 . , 1 / JO. , 21 3 J . ,270 9, , 3200 . ,3700 . , 4230, , 10*3 . 
TRXH (1, 8» =500. ,500. , 700 . , 1 000 . , 13ir5, ,16 00, t 2000 ,,2400. ,2800. , 32Q0. , 10*9. , 
T7XH (1,9»-2*;*D ,, 

I LEM =1,1, 1,1, 1,1,1, 
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SMINP 

EPSF-i.E-S*£PSG3 = l.E-5 r . 

send ; ' . , . ; ' . 

SMINPUT 

6 . 0 . 53 Z . 2 1 _H = 1 0 0 G G . * S— 17-0^* ^ w- 1 3 3 3 Q . ■ , r. =4 GOOG^-, XLE = - 1 « 
tMftX=iOQ . ,DELT=-i. , ‘ — 
QB A* -4 3 0 3 * , I C 3 = 3 , X M3 =6 .'7$n*X * B AR= 3 . , C MC L C= 1 » 5 , 1?. 
ITS- 3 » HOL-T = 1 0 . , DM A X=l-7~0-G 0 J , ♦ I MCR=53Q , EMI N = 8 0 0 0 3 .f 
W=9 5 4 $ . , X N 3 = 7 . 3 3 , C.9 A R = i 7 0 u . , X 4 B A R- i . 7 2 , 
EPSCF=.if-l-3ilN3=i,f4t>RN=5,MfQ = 6,HT=G.t 

-xaAR:=5.» : rv 

NC=30, 

DLLME=1GGG0 . , X LX = 0 E*2 0 0 J . , XLE--. 97500 37, 

SEND - * 

SMINP2 £. ■ - - " ■ 

S=53 G. * .1 “3 7523 ♦ , X = 3 • » IN CR-3 OG , CNCLC-1* , Q 3 AR=2 IS 7« 
XNB*7. »IG3=i»H^25a0a.- t ITS=3,I'iTd=l»naiT = 5D . » 

EB3- 13 fl 3 9 G • , 

-XkS=~2 . 3 9 63 7 75 Gi 4 3 1 4 7 * E = 6733 5 . 7 36 > 4 8 69 $ 7 , QM A X-2 3[ 

XLX=-.47o5632016 7P54E5,EMIM=120aGU, f 

XLX = 0, , £=2262. 727, XL.E=- .9-.U427, 

IN3=1, 

$ E NO 

Stf INPUT .... . : 

-f«Ti-3t ; . 

SEND- •• - 

sminp ^ :• ■■■■ -■ 

LT = 4 t : 

J-E-fltG 7 — - ; -- — 

SMT&&IJT , • - ' . .',V 

lNTB=t, r - .v:; 

•• . -r 

3MIMP2 - ■_■ ■ . 

SfcMO 
iMIMPUT 
I NT 3 = 0, 

SEND 

$MINP ; 

LT=G, 

send - ” ' '■ ' 



. ,XL_X=-i ? f 

!T‘1=1, £33=130339., • 
i , PSG S= • 1 e - 1 o , . . .. , 


325, XM5AR=2, 2, 
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OUTPUT, ENERGY-TURN PROGRAM 


Time history of a trajectory. 


1. Aircraft One 

IND = 1 Normal printout at print times, TMAX, and whenever ETA changes. 

Time, T , in seconds from T q ; energy ; X angle in radians ; A^ ; 

i altitude; velocity ; Q; Mach No. ; tan p; p (density); thrust; 
HM (altitude from ONEDIM ; HPRM (Hamiltonian); XDIF ( X A “ Xg ^ ; 
E; x» A^,; A^; ETA; A; PDV 

IND = 2 Additional printout for debugging purposes (see PRINT subroutine) 


2. Aircraft Two 


Same as for Aircraft One but without XDIF 


Output Notation 


When IND ^ 1 , output is 


T ~ 

E 

X 

LME 

LMX 

ALT 

V 

Q 


time in seconds from T 

o 

energy 

X in radians 

A„ 

E 

A 

X 

altitude 


velocity 
q (fcpV 2 ) 





M 


Mach number 


TANU 

RHO 

TR 

HM 

HPRM 

D '• 

l 

SINU 

COSU 

CDO 

CDB 

CLB 

CDCL2 

HR 

PHRH 

PAH 

PTH 

PTE 

EDOT 

XDOT 


= tangent u 

P v _ ■■ 

=s thrust 

« altitude from 1-D search 

» augmented Hamiltonian 

- ' - drag 
= lift 


sine jj 
cosine p 

zero lift drag coefficient 


speed brake drag coefficient 



max 



nW 

Sg 


Induced drag factor 



augmented Hamiltonian for first guess of 1-D search 
bp/bh 


ba/bh 

bTR/bh 

bTR/bE 

E 
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XLED 


A 

XLXD 

- 

• 

X 

X 

ETA 

- 

V 

A 

ss 

speed of sound 

PDV 

= 

bD/bV 

PCDOM 

= 

bC^ /bM 
o 

PCDBM 

= 

bC /bM 

d b 

PCD CM 


be /bM 

\ 

PCLBM 

= 

bC L /bM 

B(l-5) 

= 

0(1-5) 

HH1 

- 

augmented Hamiltonian from 1-D search 

HH2 

- 

augmented Hamiltonian on /3 boundary 

4,0 

HH3 


augmented Hamiltonian on /3 boundary 

i j OX* o 

EXTRA 

= 

extra term in X_ equation 

Oil 

HMZ 

= 

H ° n As, or 3 boUndary 

HMQ 

= 

H on ^ boundary 

PHU 

- 

bH/bpi 

PHH 

- 

bH/bh 

UK 


indicator for boundaries 

1 - j3 boundary 

2 = 6 boundary 

a 

3 - j3g boundary 

4. = ^ g boundary 
5 = 1-D search 



SAMPLE JOB SETUP, ENERGY-TURN PROGRAM 


Setup for CPC 66QQ 


All control card punches start in Column 1. * 

ID12345, CM75000, T300. ID No. „ amount of core to be used, 

total running time. 

PtUN(S) Set to run. 

SET(O) Sets core to 0. 

LGO/REL, LL=777777/. Will allow additional lines of output 

to be printed, over and above normal 
job cut-off. 

7 

8 End of record, 

9 

Source Deck 
7 , 

8 End of record, 

9 

Data 

6 

l End of file.' 

9 


ORIGINAL PAGE B 
OP POOR QUALITY! 
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SUMMARY OF SUBROUTINES, ENERGY-TURN PROGRAM 


MANE Main Program 

Calls PREPR 

PREPR 

Reads input. 

Preprocesses when necessary. 

Calls PROG or PROG1 
Calls EXIT 

BLOCK DATA 

Contains data common to the entire program. 

DBLSPL 

Called by PREPR 

Calculates moments to be used in spline fits for thrust and 
partials of thrust with respect to altitude and energy. 

DERIV 

Called by PROG and RKCAL 

Calculates derivatives for equations to be integrated for Aircraft One. 

DERIV 2 

Called by DERIV 

Calculates derivatives for equations to be integrated for Aircraft Two. 

ERROR 

Called from various parts of program, namely: 


FNDIX 

ID No. 40 

FNDHM - 

U J 

FNALT - 

" 100 

SLF 

" 111 or 112 

SLF2 

" 121 or 122 

SPINT 

» 334 


Writes ID No. for error location. 

Calls PRINT 
Goes to next case. 

FINDH(Hl) 

Called by DERIV 

Finds the altitude as a function of A , A , E. 

Ju* V 
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FINDH2(H2) 


Called by DERIV2 


FNALT 

- . - Called by FINDH 

■ Used when ITS - 3 . 

- Scans altitude range for minimum Hamiltonian and chooses 
proper boundary to follow. 

FNALT2 

Called by FINDH2 

FNDKM(E) 

Called by FINDH 

Finds altitudes which cause £4 5 6° to zer0 

and chooses the desired altitude. 

FNDHM2(E2) 

Called by FINDH2 

FNDI X(Bi> N, X, XTB) 

Called by FNDHM and NCALC 

Finds upper index IX for an interval in an. N point table XTB 
within which X falls. 

INPUT 

Called by PREPK 

Reads input by means of namelist for Aircraft One. 

INPUT2 

Called by INPUT 

Reads input for Aircraft Two. 

MATIN (A, N, B, M, KEY, DETERM) 

Called by DBLSPL, SNGSPL, ONED, and ONEDIM 
Solves the system Ax = B . 

A is an nxn matrix, B is an n-vector. 

Answer X covers up B , and A is altered. 

Method used ls a double-pivoting Gaussian elimination. 

If KEY 4- N , matrix is singular. 

DETERM contains determinant of A. 
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NCALC 


Called fay SLF 

Finds thrust (TR), bTR/bh, bTR/bE, A, bA/bh,CDO, bCDO/bM, 
CDCL2 , 5CDCL2/ bM , p, bp/bh, CD, bCD/bM, CL, and 
bCL/bM by using cubic splines. 

NCALC2 

Called fay SLF2 

ONED(ITT, NTAPE, CAYMIN, EX1, MPTS, 0, CCAY) 

ITT - Counter for number of function evaluations 

NTAPE - Output tape number 

CAYMIN - Answer 

EX1 - First guess 

MPTS = 4 

CCAY - Initial step size (Ax) 

Called by PROG1 

Finds x for minimum F, where F-f(x) evaluated by 
SLFX ( F, EX{ 1)) . 

This routine is used to find the X which will find minimum 

HI( = < H 2 < \> - H !( " H*) 2 ■ 

ONEDIM{ITT, NTAPE, CAMIN, EX1, MPTS, 0, CCAY) 

Called by FINDH 

Same as above, except this subroutine is used to find h for 
minimum Hamiltonian where H = f(h) . 

Calls SLFUN(F, EX(1)) instead of SLFX(F, EX(1)) . 

PRINT 

Called by PROG and ERROR 

Prints out for Aircraft One as explained in the Output section. 

PRINT2 

Called by PRINT 

Prints output for Aircraft Two. 

PROG 

Called by PREPR 

Calculates a trajectory from T ~ 0 to T “ TMAX or from 
T = TMAX to T = 0 depending upon BELT. 

DELT > 0 forward 

DELT < 0 backward 
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PR0G1 


Called by PREPR 

Sets up a table of 9 X T s, starting with X =0. 

* X 1 

Does a scan of HI( X ); whenever HI( X ) changes sign, 


X for HI{ X } = 0 is estimated by linear interpolation, and 

X X 

the estimate is used as a first guess for a ONEDIM search on 
2 

(HI( X )) . The first two crossings are recorded and, for 

X 

each solution of X obtained, a trajectory is run. 

X 

Inputs needed for these iterations that result from computations 
in the hodograph program are: 

Evader 


XLE - PAR-LFT-CEIL 
TO ENERGY 

E SPECIFIC ENERGY 
XLX CORNER MULT 
RATIO-EV 


Pursuer 

XLE - PAR-LF-CE-PUR From Loft-Ceiling 

TO ENERGY Matched Rates Table 

E PUR-EN 

MATCHED 

XLX CORNER MULT 

RATIO-PUR 


RKCAL(ID) 

Called by PROG and PREPR 

This is a fourth-order Runge-Kutta integrator with a round-off 
control. 

ID = 1 Initializes the routine. Calculates coefficients to the 
full capacity of the computer. 

ID = 2 Normal integration step. 

ID = 3 Restart - Q block is zeroed out (round-off control). 

SLF 

Called by SLFUN, FINDH, FNDHM, and FNALT 

SLF2 

Called by SLFUN, FINDH2, FNDHM2, and FNALT2 
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SLFUN (HH, RR) 

Called by ONEDIM 

Evaluates HH where HH = f (h + RR * RDL) . 
h Is in COMMON* 

RDL - +1 or -I and is also in COMMON. 

bH/bh and bH/bjj, are also calculated when RR = 0. 

As a by-product, we also get TANU, SECU, COSU, SINU, 
D, and TR . 


SLFX(F, X) 

Called by ONED 

Evaluates F = f (x) = (H2 - HI - H*} . 
X = A 

X 


SPINT 


SPLINE 


Called by PROG 
Interpolates for A^ = 0 . 

Called by NCALC 

Computes cubic spline fits for various tables. 


SPLN 

Called by NCALC 

Computes cubic spline fits for double table look-ups. 


SNGSPL 

Called by PREPR 

Calculates moments to be used in computing cubic spline fits 
later in the program. 
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TECHNICAL SUMMARY, TANDEM-MOTION PROGRAM 


The "tandem-motion" program calculates time histories of specific energies 
and the corresponding adjoint variables for 2-D motion of two aircraft at a common 
altitude which is the "loft-ceiling" (highest altitude for vertical equilibrium in level 
flight) of the one identified as the pursuer (Refs, 7, 8, 9, 10). Energy modelling and a 
thrust-along-the-path assumption are employed as in the energy-turn program. 
Throttle is full forward ( 77 = 1 ) and wings level (p,” 0) for each aircraft. 


Subscript 1 denotes evader and 2 pursuer. It should he noted that A ^0 

ilil 

and X„ 0 ^ 0 for full-throttle flight in the differential-game formulations, i. e. , the 
signs are normally opposite . This affords unusual opportunity for error in transi- 
tion between the two-aircraft energy-turn program and the tandem-motion program, 
for in the former the signs of X„„ and X„„ are the same for full throttle, both 

Jiil kiZ 

negative, this to permit use of the program for either aircraft singly in the same 

way as in the single-aircraft program (Bef. 6), The following adjoint equations 

account for the constraints h = h - h that specify both altitudes equal to the 

1 2 

pursuer's loft-ceiling. 

V 5T T g 

V „_ x ;ALi + _ii 

A E1 A EllW 1 bE^ V 1 W 1 

d CL 
^C 2 

W g L 

I s r 1 
q^ a L diV^ 


1, S,g 


lOlf I 

VC La 


T, 


ol 


dM. 


A c 

V T) 

1 ol 


] 


V 


]} 


MGE IS 
OF POOR QUALITY 
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V bT 
r 2 2 

T 2 

s q 2 s 2 s 

dC i 

02 

L'O -' | 
O' 

V 2 W 2 W 2 U 

dM 2 

d v 



W s 
2 13 r 

L 2 

2 11 


q. S n a L 

d M. 

1 

1 

1 


h 2 2 

2 

2 


. dh 2L 

A E2 d E 


6T 2 S T 2 

q 2 

f 

1 w 

bli V W 

2 E 2 2 

M 

^ 1 

2 

& 



Sj 


2 02 


^ S q 2 S 2 

V' V 2 W 2 ' 2W 2 dlJ 


d v 


L 2 ( g , M 2 da 


hW 

dM 0 

4 j 

UV 2 ' 

a dh 

, dh 2L 


6T ! 

gX, 

A E1 dEg 

l w 1 

6h l E. 

V 2 W 2 


, r JL'ni 

+ S c 2 v P dh ' v 2 y iJ 

L 2 2 

, s v dc n 

1U ( &' d a \ ol 


Si 


, 3g q x S x S x V 1 d p , 


) v V W 2 W dh i 
ol 1 1 1 


+ _1 da N /I dp _ 

V aV, a dh / U n o\p dh TT 2 7 J J 


^1 S 1 V 1 L dM l ^ aV l a dh ^ ) c 2 v pdh V J 

L 2 
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INPUT SPECIFICATIONS, TANDEM-MOTION PROGRAM 


The tandem-motion program has basically the same inputs as the two- 
vehicle energy -turn program. 


A. Inputs for Preprocessor 

namelist/minp/atab, hatb, cdot, xmot, CDCLT, XMCLT, 
RHT, HRT, TCLH, YMCL, TRXH, XMTRT, HTRT, NMA, NCDO, 
NCDCL, NRH, NCB, NXT, NHT, IDEN, IDG, LT, NEE, NEI, ETB 


IDG 

NEI 

NEE 
ETB 

All other terms are explained in the two-aircraft energy-turn section. 


Notation Explanation 

= 1 Generates power ceiling (h g ) and loft ceiling (h.^) 
tables versus energy. 

= 0 Does not generate tables as they are already available. 

Number of points to be splined in the power ceiling and 
loft-ceiling tables for Aircraft One. 

Number of points for Aircraft Two, same tables as above. 

Twenty point table of energies to be used when generating 
above tables. 


B. Main Inputs 
1. Aircraft One — - Evader 

NAM ELIS T/MINPU T/GO , S, W, E, X, TMAX, DELT, NEQ, NPRN, 
ITB, INTB, CNCLC, XNB, ICB, XLE 


■ •• .:•! 
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Notation Explanation 


ITB = 1 Uses XLE input as initial X for integration. 

E 

2. Aircraft Two — - Pursuer 


NAMELIST/MINP2/S, W, E, CNCLC, XNB 4 ICB, XLE, EDOT 
All other terms are explained in the two-vehicle energy-turn section. 

Main inputs $MINPUT and $MINP2 must be read in at least twice for 
the first trajectory. The first time is needed for the preprocessing of the h^ 
and h^ tables. 
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C. Sample Input, Tandem-Motion Program 



35, .02115, .02135, .0210, .0 210, .0242, .03241*0 359, .0 374, . 0 33 4, 
• -■ -> -• •} » ■* w j • 7 bu397 , * J 403,. 0 403, *0403,0. , 

0 1 G-»t • 2 , . do 9 : j -» -9 5-ir i. » 0 0 ,-i-» 05, S. * 1-0', 1 » 2 il,-i >40,-1 *00,1.30,2*', 

! , 2 . 4 , 2 . b , 0 . 


ShtltL— 



ORIGINAL PAGE IS 
OP POOR QUALM 
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SMI NP 
Lr = i , 

NMA=2 0 , N 9 H = 2 0 , 

A T A 3 - i 4 -} & . , 1116. , 13 4 7 • , 1077.,lul>/ , «,i0 37. ,101 E> * , ' } ' j '• » 7 ,0 38.1,9 Ei.l-j.l'ifl.l, 
963. 1, 363.1,983. 1,968. I , 968. i , 96 3, 1 , S84 . 2 , 1 Q H 4 . , 1 ;.:5o . 4 3 , 

HArd--i.E5,G.,5.E U l.E4,i.5E4, 2 . C4 , 2 , 5E4 , 3 . E 4 , 3 . 63 9E 4 , 4. - 4 , 4 . 6 - 4 , 3 . E4 , 
5.5t4,6.£4,7.E4, 8.E4 f 8. 2 02E4 , 3 . C4 , 1 . E5 , 2 . h 5 , 

P.HT = 8957 • E -6,2377* E-6, 2048 *£“Si i /75.E-o»1496. E - •>, 1 2 ab ♦ t -6 i IGA^.E-t, 
n -E33.F -7, 7061 .£-7,5351 . --7,4601 .E-7, 3613. f>7, 2 346.6-7,22 33.E- 7, 

1334 »E- 7, 85 5 6. E- 8, 7764. €-0, 51 91. £-8, ,jj. 31. £-3 ,9 4. £-3, 

M9.T — “1 * c, 5 j D * 3,1 »E4,1.S - 4,2, E 4,2 »5 l 4,3*E4, 3* 4,4 .r4, 4* 5F4 ,p. !_ 4, 

3t5r. 4)3 t z4>7 • L4,f},£4,8,2^J : ■ , 4,9*F4,i.h..y,2.z.6, 

N C 6 = 1 6 , 

TG)_4— .8'ji .•3'. , ,»3?5,»73,»76 j» 74, .8 ’-t ,1.0 4,1. 2,1*27 
4 * 0 » i 

MCD0=15, 

X M 0 T — t] • t t ? » . 1 j • 6 , * 7 7.87 »9,.95,1. ,1.05, 1*1,. 1.2, 1*4,1. 4,1*3, 3, j 4 ^ G • , 

CD3T = . 0 2,.02i.019&,.C192,.aiay,.Gi89i.02G3i.3?7,.:3^,.0:»a5i.a436,.04G3i' 
.041 4, .8413-, .G41.1, . 0413 ,4*0. , 

NCOS L= lit 

XHijLT-dt 7 j^j 4 4-), 6, , •S,.4 lT4 -3'>,ii»1.0 4,I«L,l»2»L»4»l*t,,l*8,3*,4 9 j., 

C l >6 L T — o • f}3»4jb*4io«H7 0*ftj7 >6t7 • 1 ? ■).t3,6.4,',*3*-^.6,3*8,3*,2*— ,2.6, i ^ 0 . , 

NXT=10,'4HT-9, 

X 4 T 2 T — 0 1 j . 2 , » 4 , • 6 , •81I* 7 1*2, 1.4, 1*6, 1.8, 1 U 0 » 1 

HTRr =- 15-10 0 . 7 0 . , 13 13 i3 0 * 1 2 900 0 . , 36 089 , ,40 0 01; , , 4 3:j ^ * , 5 GOG j . , 175004* , i 1 * .] . , 
TkXH ( 1, 1) = 910 0 .,9100 . , d 30'* 1 107GC* , U5Q0 .1 12101 * , 132130. , 14.530- , 1=5410* , 
iE5u a • , io* u . , 

T KX H tl , 2 ) -ooO'i • 1 6 0 u * ,7 2 02 • ,7900 * , & 7 0 1 * , 96 9 0 • , 1 ■ 1 3 ■ > u » ,11130* , 11900 • , 

12702. ,10* , 

TPXH ( 1 , 3) =4 10 J . , 410 0 . ,4 5U0 . ,51 Q 0 . ,5900 . , b ■ ? Q : J . . 7 Ml* . , 7900 . ,840 0 . ,830 0 . , 
10*0., 

T«XH a , 4) =2 9 0 U . , 2 9 0 0 . , 3 ) (JO . ,36 UU . ,4100. , 9 0 00 . , 9 '0 . * , 69 u 0 * , 78 U0 .,3203*, 

10*9 . , 

T-sX H ( i , 5) -10 0 C » » 1C 0 J . ,1 3 00. ,22 0 0 * • , 2 o 0 G * , 3300. , 4 L 0 : . * , 4 .8 0 0 , , 6 4- 'J 0 • ,5300* , 

10*0 . , 

T-<XH a ,6) = 700 . ,700. .,1203 . ,1700 . ,2700 , ,28GU. ,340 3. ,4000. ,4600 , ,5203 , , 11*0 . , 
T.9X2 Cl, 7)-dC0. ,300 . , 113 3 .,1400* , 1700. ,2100 * , 2 73 3. ,32 GO., 37 00 . , 420 u. , 13*2. , 
TSXH Cl , 3) = oT C . ,30 3 . , 7 0 0 , , 1 00 0 . , 1 3.10 . , 160 0 . , 20 3 1 * , 1400 . ,2 40 fl* , J?0G . , 10*3. , 
T4XH ( 1 , 9) =2U*U 
IDE 4=1, 1,1, 1,1, 1,1, 

IfMD 


ORIGINAL} FAGB IB 
OF POOR gUMffil 
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$hiNP 

IDG=:UN£E = i2 T LT=3,rt£L=12, , 

N££ ='13 
St NE3 

$ MIN PUT • . ■ - - 

S=50 0 , , W-36 D Q Q . , I NT3 = i , GMr.LC = l - 5 ? 

W-9545. ,XNB = 7,3o ,IC8 = 1, S = l70i , - 

&6= 3 2* lot - ; v‘: ■ — r • .7 • 

T MAX = 6!J 0 .1 NPRN -1 1 , OE L T = . r ME3= 2 » 

TMAX=90u*i 

T MAX = 3 0., - ■■ • 

SEND .- v ■■ •. ■ ■» - - ' 

SMINP2 

I.CB=4T-S=5^a,»M=37523--,XH3=7i, 

CfiCL G- L . , 

SEND 
SMI H PUT 
£=45000. i 
SEND 
SHI tip 2 
E=46Q8Q .35* 

SEND 
Smew put 

I NT 3=0 t 
SEND ■ 

SMI HP 
LT=0 j 

scyo 
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OUTPUT, TANDEM-MOTION PROGRAM 


Evader 


T 

E 

PHE1 

ALT 

VEL 

Q 

MACH 

MU 

RHO 

TR 

D 

L 

CDO 

CDB 

CLB 

CLH 

CL 

CDCL2 

EDOT 

HDIF 

A 

BET7 

OME 

H7 

H6 

NEQ = 4 . 

the following 

is printed. 



LME1 

LME2 

CAPLH 

LME1D 

LME2D 


Pursuer 


T 

E 

PHE2 

ALT 

VEL 

Q 

MACH 

PHEE 

RHO 

TR 

D 

L 

CDO 

CDB 

CLB 

CLH 

CL 

CDCL2 

EDOT 


A 

BET7 

OMP 

H7 

HG 


PHEE is analytical (closed-form) { bh /bE ) printed as a check. 

i £t 
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ABSTRACT 


This manual is a guide to the use of the AIRCRAFT ENERGY- TURN HODO- 
GRAPH PROGRAM which produces an extended energy-maneuverability description 
of the capabilities of an aircraft whose characteristics have been input. The 
various quantities computed are those needed in the energy-turn and differential- 
tur"' analyses of the references. The program is in FORTRAN IV for use on a 
CDi.': 6600 computer but with minor modifications can be ran on any computer ac- 
cer , ng FORTRAN IV. 



INTRODUCTION 


The AIRCRAFT ENERGY-TURN HODOGRAPH PROGRAM is a computer 
program for calculation of aircraft maneuvering performance in terms of the 
following quantities as functions of specific energy level: maximum energy rate, 
maximum instantaneous turn rate, maximum sustainable turn rate, power-ceiling, 
loft-ceiling, energy rate at loft-ceiling, and loft-ceiling rate at loft-ceiling. A 
choice of approximations for treatment of vertical thrust component is available 
on input option. 

The program is related to the SUPERSONIC AIRCRAFT ENERGY-TURN 
COMPUTER PROGRAM of Refs. 1 and 2 and requires the same aircraft tabular 
data represented via cubic-spline lattices (Ref. 3). The coding takes advantage 
of the NAMELIST feature of FORTRAN IV, which permits 'stacking 1 of succes- 
sive cases to be computed with minimal input. The program may be used in 
conjunction with the two programs of Ref, 7, as described in Ref. 8, to generate 
families of differential -turn trajectories and their barrier surfaces. 
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ANALYTICAL DESCRIPTION 


Turning duels between two dissimilar aircraft are studied In Refs. 4, 5, 
and 6 with the use of energy modelling approximation and a differential game for- 
mulation. Computation of various energy-rate and turn-rate quantities are needed 
for comparisons which represent steps in the numerical solution of the differential- 
turning game, and which are in themselves of more than qualitative interest for 
comparison of configurations. 

Power-Ceiling at a specified energy level is defined as the highest altitude permitting 


horizontal-force equilibrium in level flight. Designated h , it is determined im- 

o 

plicitly by 

' 2, 

0 

with 0! from 

(la) 

qSC a + Ta = W 

Li 

a 

o 

(2a) 

The computation is performed with Specific Energy , E = h + V^/Bg 
An alternate thrust-along-the-path option is available with use of 

held constant. 

T - D = 0 

(lb) 

qS a = W 

(2b) 


CL 

Either eq. (la) or (lb) is solved computationally by Newton-Raphson iteration be- 
ginning with a first guess determined by a rough scan over the altitude range. A 

quadratic C versus C model is employed, with no upper bound imposed upon 
D L 



Loft-Ceiling , , is defined as the highest altitude permitting vertical equilibrium. 

It is determined as the root of 


q S CJM) + T sin a* - W (3) 

is , a function of Mach number. Stall angle -of- attack, a * , 

max 

may be input as zero for thrust-along-the-path modelling; this choice also triggers 
the use of (la) and (2a) in preference to (lb) and (2b). 


where 


C L (M) 


Sustainable Turn Rate is computed in small- a approximation as 


io 

P 


g sin u 
VW 


(L + Ta) 


(4a) 


A candidate a is computed from the horizontal-foree-equilibrium eq. (lb) ; 

employed if the corresponding lift coefficient 0^=0^ a £ C^ ; otherwise 

a 


it is 
C L 


and corresponding a are used in (4a). Computations are performed at fixed E 
and h . The maximum sustainable turn rate for each E is picked out by a scan 
over the altitude range bounded below by terrain limit, , dynamic pressure 
limit, h f or Mach number limit, h , and bounded above by h . The altitude 

Z tJ O 

bounds h and h are determined by iterative solution of the defining equations 

u u 

jS = 0 and 0 L- 0 as in Refs. 1 and 2; in fact, the coding employed is identical. 

Z o 

The approximation to the maximum sustainable turn rate found by the scan over a 
coarse mesh is further refined by the use of golden -sect ion one-dimensional search 
starting at the altitude mesh-point chosen by the scan. 


An alternative estimate of sustainable turn rate employing a thrust-along- 
the-path model is triggered by input of a* = 0 . (The estimate of the preceding 
paragraph is automatic for a* f 0 , as is the power-ceiling estimate of eqs. (la) 
and (2a).) It is important to input a* - 0 when compatibility with the program of 
Refs. 1 and 2 is needed. The thrust-along-the-path model uses cos a = 1 , 
sin a = 0 , obtains a candidate as 
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( 5 ) 



and employs it if < CL , otherwise C in 

J-j L 


W = 
P 


g sin ji q S 


VW 


(4b) 


Maximum Instantaneous Turn Hate is computed as that for at the corner- 

velocity altitude, li , if h^ 5 s o for the chosen energy level, otherwise at 
h = 0 . The corner-velocity altitude is determined by iteration on 0 - 0 and 
jSg = 0 (C^ limit and normal-load-factor limit) just as in the program of Refs. 1 
and 2. 



£ 

V 



+ T sin a* 
W 


2 

N 

) 


1 


( 6 ) 


Energy Hate is computed in small- a approximation as 

2 

E - V [ T ( 1 "f“)- D ] ( 7a ) 

With a determined from eq. (2b) at each specified energy-altitude point. This 
approximation is used when a* ^ 0 is input whereas, for c£* = 0, a thrust- 
along-the-path model is used for which the corresponding equations are 


E = V[T - D] 

and a determined by 


CL 


W 


(7b) 


(2c) 
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Standard print-out for both options consists of the maximum energy-rate value at each 
specific energy level, determined by a scan plus search between bounds identical to that 
for maximum sustainable turn rate, and the energy rate at , the loft-ceiling. Also 
printed is time rate of change of loft-ceiling altitude. 

TURN RATE COMPARISONS 


When the characteristics of two aircraft are input, computational comparison of 
sustainable turn-rates is carried out. One such computation determines the curve in 
energy space separating regions of sustainable-turn-rate superiority of each craft over 
the other; the other compares sustainable turn rates at matched loft-ceilings (Refs. 4 
and 5). Prior to these computations, cubic spline fits are made of maximum sustainable 
turn rate versus specific energy and loft-ceiling versus specific energy for both aircraft. 
The computations then proceed iteratively for several input values of sustainable turn 
rate, yielding as many as four pairs of energies for each. It should be noted that the 
equal-sustainable-turr.-rate curves in energy space enter the determination of the number 
and sequence of subarcs of trajectories and affeet the barrier surface structure. 


A similar computation of pairs of energies for equal instantaneous turn rates is 
available on option. 


TANDEM /LOFT MOTION APPROXIMATION 


A rough estimate of the tandem/loft motion starting at the matched-loft-ceiling 
curve and proceeding for a single step backwards in time is provided by Euler extrapo- 
lation. If (h - h ) £ 0 at loft-ceiling' match, tandem/loft motion is a candidate option 

for the pursuer and, accordingly, the computation is performed; otherwise it is by -passed. 

•« •» 

If the second derivative of the rate difference (h - h ) is positive, the time step 

J-Jii J— 1 1 


At 


^ 2 ^ 


( 8 ) 
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is calculated and employed so long as it does not exceed in magnitude an input value 
At<0, chosen as perhaps -20 sec. or -30 sec. If the second derivative is not positive 
or if At given by (8) exceeds At in magnitude, At is used. The pair of extrapolated 
energies and the At used are calculated and printed for various initial points along the 
matched-loft-eeiling curve. 

CHATTERING SINGULAR ARC 


The locus of a singular solution in energy space (Ref. 6) is determined upon 
input option by iteration on A = 0 , A being a determinant having pursuer and 
evader energies as arguments. 


A - ( £ 2L - f 2)( f i f Tl 


12L 


M bf 

-f + 

bE ML bE 2 bE J 

i X & 


" ( f 12L " f l ) ( f 2 


f 

2 bE 2 - 2L 


Of 

-A) 

bE„ / 


= 0 


where 


2L 


rWLb 

i 

f* = 

V (T - D ) 

r 2 2 1 

L~ W 1 J 

h i* 

w J 
2 1 

rW I V-| 


- = ! 

r T i<V D i>-| 

L w 2 J 

h„ T 

’ X 12L " l 

L w 2 J 


2L 


O' 

hC 5 

t 

b 

— i 

3 

1 

d 

H 1 

I 

1 + b I 

-V^-Ei) 

dh^ 

bE i 

bE 

L . 

L* Bh 1 
h l 1 

J 

h ; dE i 

bf 2 

b 

- V 2 (T 2' D 2 ) -] 

b 

l . 

rWV 

* 

i dh 2 

bE 2 " 

bE 2 

- W 2 J 

s “■ 

L w 2 

h 2 ^ 

bf 12L 

b 

r v i' T i - D i 

-1 



bE_ 

1 

bE 

1 L W 1 

J. 

h 2L 




bf 


12L 


bE, 


a r v i< T i- D i> i 
bh^ w J dE 2 

2L 


2L 
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bf 

2L 

b 

- V 2<V D 2>- 

1 _A_ 

- V 2<V D 2>-| 

6E 2 

- bE 2 i 

' W 2 ' 

J h 

h 2L 2 

- W 2 J 


bh 


_2L 

bE„ 


‘2L 


In these expressions, h. is the altitude for maximum energy rate f. , and h, 
is the pursuer’s loft-ceiling altitude. 


2L 


The locus of points for which A vanishes subject to the evader's loft- 

ceiling exceeding the pursuer's is of interest. To determine this approximately, 

a scan is carried out of A versus evader energy at fixed pursuer energy, using 

the loft-ceiling match energy previously determined as a lower bound, A more 

2 

precise determination is then performed by minimization of A via golden-section 
search. This is done for several input values of pursuer energy. 

Since the computer program employs cubic-spline fits of h '(E) and first 
derivatives of these spline representations, care must be taken when h*(E) for 
either or both aircraft exhibits a jump discontinuity to carry oi t the chattering 
arc computations sectionally, i. e, , not requiring the spline to approximate the 
discontinuous function, so as to avoid undue numerical error. 
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INPUT SPECIFICATIONS 


A Input Tables 


NAM ELIST/MINP/ A TAB , HATB, CDOT, XMOT, CDCLT, XMCLT, RHT, 
HRT, TCLIi, YMCL, TRXH, XMTRT, HTRT, NMA, NCDO, 
NCDCL, NRH, NCB, NXT, NHT, IDEN, LT, EMIN, EMID, 
EMAX, DEMID, DEMAX, ETA, IDGO , IDG1 , IDG2, IDSUS, 
EVL, ISEP, DLTI, IMXI 

IDEN is dimensioned 6 and each slot is associated with a table IP. The 
tables are as follows: 


NUMBER 


ID 

TABLE 

DIMENSION 

OF POINTS 

DESCRIPTION 

1 

ATAB 

20 

NMA 

a (speed of sound) 

1 

HATB 

20 

NMA 

h (altitude) 

2 

CDOT 

20 

NCDO 

°D 

2 

XMOT 

20 

NCDO 

U 

Mach Number 

3 

CDCLT 

20 

NCDCL 

% , 





T 

3 

XMCLT 

20 

NCDCL 

Mach Number 

4 

RHT 

20 

NRH 

p (density) 

4 

HRT 

20 

NRH 

h (altitude) 

5 

TCLII 

20 

NCB 

ft 

C L 

5 

YMCL 

20 

NCB 

Mach Number 

6 

TRXH 

20,20 

NXT, NHT 

TR (thrust) 

6 

XMTRT 

20 

NXT 

Mach Number 

6 

HTRT 

20 

NHT 

h (altitude) 


EVL 

100 


Energy to be processed 
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Other notation explanation: 


LT 

EMIN 

EMID 

EMAX 

DEMID 

DEMAX 

ETA 

IDGO 


IDG I 

IDG2 

IDSUS 

ISEP 

DLTI 

IMXI 


= 0 indicates end of run (termination trigger) 

= 1 indicates evader processed 

- 2 indicates pursuer processed 

minimum E 
internal cutoff E 
maximum E 

AE from EMIN to EMID; if DEMID = 0 , program expects EVL table 
AE from EMID to EMAX 
constant multiplier of thrust tables 

= 1 processes evader when two-aircraffc study is desired 
= 2 pursuer processed 

= 3 one aircraft study to be followed by study of another aircraft 

=4 one aircraft study to be followed with another study using the 
same aircraft 

=1 normal print 

= 2 skips table prints other than matchings 

= 1 no chattering 

- 2 includes chattering 

number of points in grid of matching maximum sustainable turn rates 

= 0 no approximate tandem loft motion 
f 0 calculates and prints approximate tandem-loft motion 

limit on A in tandem-loft motion; must be positive 

number of points in grid of maximum instantaneous matchings 
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B 


Other Input 


NAMELIST/MINPUT/GO, S, W, CNCLC, XNB, QBAR, XMBAR, EF,3, 
HT, ICB S ALPN, ALPD 


Notation Explanation 


GO 

gravitational constant (32. 16) 

S 

wing planform reference area 

W 

weight of aircraft 

CNCLC 

speed brake factor, where CNCLC 

N (maximum structural load factor) 

XNB 

QBAR 

Q boundary for (3 calculation 

a 

XMBAR 

Mach boundary for /3 calculation 

u 

EB3 

If E < EB3 , by-pass A calculation 

1 5 

HT 

maximum altitude allowed 

ICB 

= °. 6,= ! 

= 1, = TCLH (MACH) 

ALPN 

numerator for a ' 

ALPD 

denominator for 


OUTPUT 

Sample output follows. 
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rfPDiV&R.VH 9 ATE ' ATH 


»Di{£ 9 P'OW t 2- Zi I L I H '. 
CEIL I ^ ru-iN ?M- 


on5P\'f'*ri :r Al.ALV TICAL 1 r emi3'> AlSOCTATES erg. 

ALPHA'S 0 . C PURSUE •*= E4 


LCF T P*MF t - CE I L AL T J T JOE 

"£ I L I “4to tt r *j-pGY L0*rp IQ USD 


LaatP a «0 PAX SUSTNHLt fU* SUSIN 8 LE 
TURN RATE TURN RATE A.T TUVM RATE 


SPECIFIC 

ENrFJJV 


4. 71 2* WO 81 
9. 60 54 31 F +0 3 
1.44 5414 F+OA 
1.92 3ft3Jfc>0 + 
2. 3)2983 £*04 
?.84*F»19F < 0‘t 
3.285826F+14 

3. 697383004 

4. 07 17650 34 

4 . 417703 E+Q 4 
4 . 72*5 m*a<» 
4. 99654 IF+04 
5.2V)479F*C4 
5.409070E+04 
5.5SJ17*F*84 
5.61 3 4*0 F+ 04 
5.3?*999F*04 


3. 

r. 

n. 

3 .. 

3. 

a. 

3. 

J • 

0 . 

3. 

A. 

0 . 

U?2» *7SE-0 > 
9 . 46 & 416 E -01 
0 . 

0 . 


I. 751 :5?»: *1 ! 
8 . 5 * V 0 *?F» 3 J 
1.311 14 If -14 

1, '] ’-. 11 E + 14 

2 , 

2. 774*13:* 14 
3.1 7‘*2G4r + 04 

3 . blli..*& 7 *:*a 4 

4. flZ9469£+04 
4.4273.3IF*n4 

4, 79-i jo?c*04 

5 . n 4 a 3 7 r *04 
5.759526 t >04 
6. 484092E+04 
7.129955E+I14 

7. 2 7950 A £>04 
7.252192004 


9.5374? ’"-11 
9.55* »: r- U 
9.4H? 1 *9' r -i!l 
9, 371 »«*5E><u 
9.244 IhlF-fll 
9. 97*7311-01 
s. 95*4i>.r-ni 
8 . 544 <U 9t> 0 1 
8.126 vDUF -Ql 

7.734I27--01 

6. 7691 ivC-Ol 
5,312 )21£-()1 
7.7621 0 It- 0 1 
8. 258* 81 £-01 
4. 5 25563 F- 01 
i.6&&?2fE+00 
-1.34420 >£♦ 01 


0 . 

a. 

а. 
o. 
o. 

ft . 9ft 0 758'E * 32 
J.S 289 1 OF *13 

б. 0J'60 23F*03 
A.4321VSE+03 
1.0723 94E+04 
1.290 197094 
1.49 ft'H54E ♦14 
1.88VI98EMJ4 
2.234U 72F*34 
2.552477E+84 

2.*41490E^04 
3.1126956+04 


2.41436FE-01 
2. 6 534? BE -01 
2.434 790E- 1 1 
1.531275E-91 
6. 99 462 ?£ -02 
- 1 . 2 G 0 30 OE t 33 
-1. <1000 COE *33 
-1.0 GOQOOE* 33 
- 1 . OGOOOOE* 33 
-1.0 OQOQ DC *3 3 
-1.0 00 OOffiE f 33 
-1.0 00* COE* 33 
0 , 

-l.OOOQiOOE+33 
-1.Q80000E+33 
-l.OGOOOOE*33 
-1.0 000 DOE +33 


n . 
r. 

o. 

4.4629592*01 
S. 60 41 68 E +03 
i .4854 00E +04 
•’.02284 IE + 04 
2.564Z43E+04 
I .751015E+04 
7.026S93E+04 
2.283783E+04 
2.47277&F +04 
? ■ 74 59 30 E * 0 4 
3.20668AE+04 
1.772413E+04 
6.152250E+B4 
4. 5743826*04 


2.414345E-81 
2. 6 5942* E- 01 
2.434790E-01 
2.I725J9E-01 
1. W8 726E-01 
1 . 646 1 92E-01 
1. 384 94 IE -01 
1.137135F-01 
9. 85&445E- 0? 
8. 753 915E-02 
7 . 94 37S2E- 02 

7. 7216Q4F-QZ 
5.932275E-02 
4 • 748956E-0? 
3.S62927E-02 
Z.397594E-02 
9.4B1281E-13 


5.000 OOOf +13 

1.000 Q80C+Q4 
1.50 0 OOOf +94 
2.80. 300 e +04 
2.5 00 OOOF +04 

3.*0°C3-Ji *94 
3.500 0*ir+ 34 

4.000 000 r + 34 
4.500 QOOL+04 

5.00 0 OOOf *04 
5.5000'00 c + a4 
6.00j000r*04 
7. 0*0*0 or + 04 
O.OOGOOOE+04 
9*000000^+04 
1.0 00 000; +05 
1.1 03 OOOf *85 


ALT < 

r«3t-r?*!-iT 


il * 

i.497&iar*03 
5.a?58J2F*a i 
b.i42524F*oj 
i,um?t:*u4 

l.'ilSSWftBJi 
1.1181331 *04 
2.I2T 

2. 7116141 *C4 
3.0517017*04 
3.252011*' *34 
3.22804 7 r *04 
3* x** ra 7ZI- n)4 
3. 18 7977:: *34 
3.5637?4F*16 
3.666ft54f + (]4 
3.71111 316*04 


'i'AK-1 FIT. 
TUP»'-?A'' 


7. 4 ! 4 Ibbr — u 1 
v.qiz j c 3.:-oi 

2. 7<M)I5<‘-at 
?. 5*1 33F* -01 
7. 3 *04* 1C - It 

?.zvmu.-ai 

2 . 145737* -<n 

2.Q?7.»b2E-ni 
i.990I79£-oi 
1. R j >, ^70C“ 01 
I.*vi8b7l£-51 
1.4?46*rtF-91 
l.?43b53F-0l 
1. 1H22 36-01 
I* i 0 Vlfe IE * fl't 
l. 078338E-1I1 


'..1 -Oi-*Vll 


J. 

0 . 

C. 

3. t>Ti45ir*in 
9.«H>i.41t>0 T 
1.4JI j19iS*04 
t.9*J/89r >84 
1,4 b* 1 44 Jfc*94 
1*7 >2 , 'lNf>04 
?.(J1 » »80fc*94 

2. 333 LS 3r.*<J4 
2. 5-70*8 If -14 
Z . tv 1 * 1 bf»5E* 84 
3. 380 4896*06 

3. 9 a ia0.!€*04 
4 • 2 lb li 7f ♦14 
4.6lO|3ftE+Q4 


£*)■' *Vv.q\Tt 


4. 2-41 -I <12 

6.8 37434- O? 
7.25444,“# <12 
blab/^i* f *oz 
b.?*t **'K*U? 
4.4 ?0 8» ' - *92 

4.4 It H»7r * <12 

4 . 3 8 10 4 4 F> 02 
4.!24‘i4'« l .#B2 

3.9 71-140- +02 
I.914LO:r»Oi; 
4.0121? *1*32 
3, 89244 <->02 
3.57914 »T * 02 
2.9J4<* <*F* 32 
i.81S‘iij > f’ * 02 
3.9Bi5:&~F*01 


CtlL-^ATE 


4. 114*42? >12 

5. 77 794:6E+ 02 
b. 8 312 57 £*<12 
b. 435? 76F*9? 
5. 7/92486*92 
4.989?1J£*32 
4.1467577*32 
3.7435?3F*12 
3.35?l80fc*02 
3.074897E*n? 
2. 70 32 9 36 *07 
2.131 b76£*02 
3. 0213556*92 
2.95271 7E+ 02 
1. 32889 7F* 82 
3.0319726*02 
-5. 3255366* -12 


£W'= '&Y-94TE 
M-lOM-GETL 


l .55 185 IE *02 
1.45 154 1E+02 
1 .2 5 00346*82 
1.0557476*02 
8.61U20E+01 
E>.8bL3?3E*Gl 
5.0 722b 3£*Gt 
3. 16922 5E *01 
i.J;6b969E*01 
-2.644C10E*0'0 
-l.bS 1594E* 01 
-2.3i0ia4C£+01 
-5.43591 7E*01 
-i*4t?779E+fl2 
-2. 315854E+02 
-Z./7i 3666+02 
-2.9801 34E *02 


LOFT-CElL— <T 
AT-L0FT-CETL 


1 .559970E*f)2 
1 .364C58E+02 
l.lb3l06E*02 
9.60 793 6E *81 
7 • 73247 If *01 
b. 87490 86+ 81 
4.477«50E*01 
Z ,722294E*0l 
1 .11 11356*01 
-l.942485£*00 
-9.23I872E+00 
-9.1374686*00 
-5.4458056*81 
-1 .516649E+02 
-1 . 1396546+02 
2. 17 3291 E* 01 
-3.270466E+01 


S^ECIFIG-fHF* 


5.00000(11*0 3 
i.5<oaoon£*04 
1.5000 OOF* 04 
2. QCOOQOE+04 
2.5000806*04 
3.9G0OgCF»Q4 
1. 50 0008 E* 04 
4.000000E+04 
4. 500000 F* 04 
5. 000093 E# 04 
5. 500 080E* 04 
6.0000006*04 

r.aoflooact 04 

8. 00 BOO OE* 04 
9. 0 00 QOOF* 04 
1.00100006*85 

t.iaoegnt * 0 5 



■lOOOGPAfM rate data 

PROGRAMME 1 r ANALYTICAL 9EG4ANTGS ASSOCIATES TNC. 

ALPHA*- 1.0 t YAOEP= HYPO 

=»f»WE3 POAiF^-JULlN ■ LOFT PAR-CFI -C'lL ALTITUDE LOHER RNO MAX SUSTNRlE MAX SUSTNPLF SPECIFIC 

CEILING TURN R4TT IEIlING T9 T'l^&Y LOW t 3 SOUND TURN RATE TURN RATE ALT TURN RATE ENFRGY 


4.142389O0S 1. 1 . 341‘.B r -+0,j R.&8 J’?tF-0 l 0. 1.42J&5JE-01 B. 1.4 7965 JE-C L 5.09f>OOQE*93 

8.858rrOF*'u ij. 8.77> ■95**1] I 9.AJ27? ’F-01 0. 1.507697E-01 0, 1.5BT697E-01 1. 0OOO8iJF*O4 

1.346402* *04 2.2?5164t-n e 1.367i?5E*l4 9. 566941 L- 11 0. 1.5519r?E-<)l 0. 1.553972E-01 1.5OO0O9E*O4 

1.79219»F*94 1. 4?fl31tE**0*> 1 . 1 3’ 564F*14 9*4 579-3 JE- (U 2.9>0m£E» 13 1.427962E-01 ?.901O42F*3J U427962E-0I 2.0900Bf£*44 

2 .2 £ O'l S4:t. * 3:4 7. 54:9956.= -3 7 ? . 30 J/*2E*!J4 9.152^7^-11 6.136771F+13 1.281Z70E- J.1 6 • 136771 E *0.3 l. 7B127QE-Q1 2.5Ofl03OE*O4 

2.6285bOF*0:4 1, 2„7?1497E*14 9,.22477'E-n 9.23*71&£*03 1.14J428E-01 9.238706E+01 1. 143420E-01 3.000 00 DC *34 

3'* 01 9 8 33 £*04 0. S. 227ll2E*04 9. tflSL'l 1^-01 1 .22O'5«*9F*04 1.02>868E-»l 1.220589E*04 1. H26668E-01 3.500 0806*84 

3.392871E+04 0. 3.< - J646E*04 8-.690 72r£-0l 1.50 il87E*04 9.393312E-02 1.503187E*Q4 9. 39331ZE-82 4.808 OOOE*«4 

3.730018E+94 ft, 4.0384 Tot +14 5.3691452-11 1 . 771 326 E *9 4 8.78447&E-02 1.771326E+04 5.7B4476E-02 4«580000E*04 

4.07204SC304 0. 4.51 ; 3756+04 S. 1u5766,E-01 2.D26271E+14 8.293742E-02 2.O26271E+04 5.299742E-02 5.0800005*04 

4.2760 32E+ 04 0. 4.9'J8.J2bE+04 7,71047 76-11 2.268923EM4 7.51354BE-02 2.265923E+04 7.A33S48E-02 5.5580OSE*a4 

6.493859F+04 O. 5.2a*935E+04 7. 453131, E-fll 2.498394E+T4 7.420093E-02 ? .498394E* 34 7.420893E-02 6.000 OOOEM4 

4.846, 761F+0* 0. 5.999757F+G4 6. 7555?., E- 01 2.921049E+04 7.001040E-02 2.921049E+04 7.001840E-02 7.0®0800E*04 

$.246<J70F*04 0 . 6.539305E+04 6 . 0461 7&E-0 1 3.310B98E+O4 6.635 J87E-02 3.31089*E*04 6.638387E-02 8.3000806+04 

5*69 0'832E*04 0. 7.21GU 32E+04 5. 367789E- qi 3.640S15r+H4 6.127722E-02 3.640536F+04 6.127722C-07 9.0000005*04 




6L1-F3T-K VX 
[ I j T - T - *i - - T 


J.flll655*'*U 
7.ga8227 r * 31 
t.U4+Q?F*n4 

I. 487975; *J4 
i.e>*ui r Ki4 
7.7191^4 14 
<?.54!»08 &f*34 
2.84 921 36*04 
3. 1 7 #9776*04 

J. 46977lf*Q4 
3. 72 7244* *04 
4. 17 35 136* 04 
4. 57531 IF *04 
4.332 783 r * 04 


14 a " 1 I it. 

rj-'j- ■’■a. 


?.956n!9"-‘»! 
1.2711 17; -nt 

1, U .' 72 ‘j)9r - 3 1 
“*m H 4*».144r -fit 

2. A 766 Ilf -31 
**.5J 7.19« .-HI 

i71*9v>-0' 

2.2 1 47 1 1 1 . *01 
2. m5 371 -.It 
1.9 a 542l*-Dt 
l. 8 J . 3917.*- 31 
t.7*6l41C-5t 
I.6U453F-01 
l . 455 )63_-<U 

1.3 J5192F-J1 


itr-CiF- i&X 
-*t- 


L * 

0 \ 

C . 

2.101 :42t+ 11 
6 , i 77 1 -mi 
'3. Z 13 704* >13 
1.2?'Z!*33f *04 
1.5U.Ui7r *04 
i .771.1266 *04 
?,i - !2u2/If *0 4 

2. *(.. *F4 

2 .4383 34 tf *34 
2. 3?i"4 Y£ *94 

3. Il J39!5fr*n4 
3.446 534£*04 


14^ L1UH 
6\* ;;v-f?iTE 


7,918 |3 * r * 07 
4.546 , *oVs;*82 

5.9 J >47 ^r+o? 

5. OH If lF*f!2 

5. 4 .4 * Jt T* 02 

*.042"77F*0? 

4. 6 . r * 0 2 

4. 5-'tfl7‘’t:*02 
4. 554 f fl ?T * 0 2 
4.6/U 136' * 02 
4. 766 lli '**Dt 
4* #91.91 2* * 0.2 
5, 5491 5?* *Q? 

6.4 78 i 0 4F * 02 

6. 9 31 '*2 1 T* tl 2 


i.K-UJ^T tNiWY-RAIt lOFT-GElL-U 
GVrt-’4TF AT-L3FT-CFU AT-L'IFI-GElL. 


?.fl?9.l43£*T? 2, 7344 3 IF *01 2.6986376*01 
4. 379277* *12 9.36U3?E*00 9.4984686*00 
5 . 664 7 73 E *J? -1.00 *51 96* 0 1 -9. 60374 IE * 00 
5.5231 1 3 r ♦ 1 ? - J.13104«E*ttl -2.9644736+Qt 
6. 1 1 6' 19C + 0? -5.29VI46E+31 -4.9687116*01 
4 » 65 14 1 2£ * 02 -7,3839436*01 -6.8161326*01 
4. 27 70 7»6 *02 - 3,5037426*31 -8 .65594GE*0t 
3 * 9 J'5 3 ? OE * 32 -1.203750E*U2 -1.0461456*02 
3.813*>75E**12 -1.51091)46*02 -1.2812216*02 
3. 78 44 sir + 3? -1. 8671626*02 -1,5134766*02 
3.7081771*02 -2.1627236*02 -1 . 68270 2£* 32 
3.6460126*02 -2.4H522E+02 -1 .812295E*02 
3.816.1.106*32 - 2 . 9 66 7 b 9£ * 0 2 -2. 003829£*02 
3.9170206*3? -3.5652066*02 -2.1555866*02 
3.7529556*0? -4.24#676C*Q2 -2. 280600E*0? 


SPfr.TflG-EHF? 


5.0000006*03 

i.iieooor.E«&4 

1.5 00 a ODE* 04 
2.0000006*04 
2.5000006*04 
3.0000006*04 
3. 6 000006*06 
4.0000006*04 
4.500IUOE*06 
5.0000006*04 
5.5000006*04 
6. 000 JO 06*04 
7.0000006+04 
8.0000006*04 
3.0000006*04 



lull. DATA 

L p r r - r i;l T v i ; -*! A T S HE □ a A T - S 



EVAl* *-r 


PU 5 -"'' Ml X- f vliT-T^M MAX-E''!ST -T 9-h 
MATCHtM AT&- r V 


+A9&IN LF-GEIL-RT-AT L* -CEIL-RT-AT 
LFT-eiXL-tV LHT-GEII.-PUR 


MARGIN 


5. Q00030E+33 
1.0000005*34 
1.53 0'S 30 F* 0'4 

2.00 00 00*: *04 

?.5o oo oof* 04 
3. do oo oo ft aw 
3*50 0609004 
A. 0000000-04 
4, 5 9 00 ODE *04 

5. 00 0 000:004 
5,53 OOOOE*04 
6.00 00 00 E* 04 

7. 000000 £-*-64 

6. t)9 9Q0DE*04 
9.0000 aftr+04 


5. 1 )/,J 0 V ♦ 0 ) 
1 » O '? 2 2 2ir * .) * 
1.5’6?n?t0'i 

2.5 {9441E + 04 

3. 04 A J 70" + 04 
3.55il64t£*04 
4. 06965 OF *04 

4. 585 t o IE * 94 
5.in8J16E*94 
6,674493E*04 
ft. 2993 16E * 04 

7. 350566E*04 
5.205>97F*04 
9. 249777E *04 


? . -j 5fcol 5E-0 1 
J.2’1 107E-01 
3..t27 299 r -9l 
, .'»4M44F-C1 
2.676 r41t>0i 
2 . 51 799HF -0 i 
2. v7r535t-01 
7.Z347I3E-41 

z.xoa .m:-oi 

1.995421E-31 

l.899U7t-ai 

1.7SP 141 E -91 
1.501 j5 JE-D1 
1.455 063£-01 
1 .3351)75-01 


2.4549:! -tE-fl t 
3.3*13^7—01 

2 . 7 TH il T -- 0 1 
Z.5 J75-4 r ~0l 
2.3 70 71 1 i -01 

2.2 7?j*:L-0l 
Z. 291 J /“r - f) 1 
2. 11644 77-11 
2.0V4^>91 
l.966’2 r v*31 

1.7 59 ?7 65-01 
1*5739717-11 1 
1. 3 73247 f-- 91 

1.2 70 (3*5 7 - 11 

1.1 69 **,7- -C 1 


- 4 . 92 71 4:82-92 
-2.1745TSE-92 

- 2 . 51 2 tl 2 0 f. - 9 2 
-3.0'8 5b9ZE-0? 
-3. 0614132.-92 

- 2 • 45.7 1 6.1 E -9 2 
-1.69&168E-52 
-9. 7S 2 5441 -03 
-4.4 11 V81E-93 
-2.91 34335 -53 
-1.00 5 465 £-02 
-2. 12 22 002-02 
-2. 25 4111 E -9 2 
-1.941 253E-02 
-t.66‘3755r-52 


2.698937E+01 
9.49A468e*00 
-9.60 3243E*0fl 
-2.9644 7 3E *91 
-4. 95371 IE* 01 
*6. 8 161 322*01 
-3,655940E*0t 
-l.il 46145E*02 
-1.29 122 IE *02 
-1.51 3476E+Q2 
- 1.6627 02E* 02 
-1.8 12295E+02 
-2.00 3629E*D2 
-2 • 1 55586E *02 
-Z. 7 8 06 COE *02 


1.553106E+02 
1 .355119E*02 
1 .1523022+02 
9.479249E+01 
7. 59576 6E *31 
5.921 879€*01 
4.28O563E*0i 
2.4779012*01 
8.6617282*00 
-4.3190722*00 
-6.956855E*00 
-1 •3970982*01 
-8 • 97432 62 * 01 
-1 .64679O£*02 
-6.157363E*0l 


1. 2 8 13 032* 02 
l. 260 134E* 02 
1. 2483 JS£*02 
1. 244 J72E+ 0 ? 
1.2554562*02 
1. ?73601E*02 
1. 293650E* 02 
1. 2939 J5£*02 
1. 36803BE+ 02 
1 . 4702S5E* 02 
l* 59 J113E+ 02 
1.672585E*02 
1. 106397E*82 
5. 3.879562*01 
1.664 66 3h* 02 



mail, in : ifiR. i’H ii A f 1 
l. I -*7r TL T .-MAT Of r> rfil- 


tVfljrs.f 


pu <-»**•< Lj-r-ef.iL sus-m-M-tv 5 us-tm-9T-pu7 

tatGHC i 


MARGIN 

MAX-LOFT- 

HAX-LOFT- 

HARCfN 


CEIL-RT-EV 

GfcIL-^l*PUR 



s. ooQOOor+oj 
i.oo3ooje+34 
1 .5000001 + 0'4 
2.39 DO 00 e* 14 
2.50 OOOOF+04 

3.ao oooq,f+q4 
3. 50 GO OOP +04 

4. 00 GO OOP +0 4 
4.50QOCJE+94 
5.0000 OOF+04 
5.50 00QQE+04 
fa.aoaooGE+04 
7.003000C+04 

5. 00 000 0 £ +04 
9.000BOOE+14 


5. 1H7J00r +Qi 
1.0’2J29 e.+ 'J4 
1 . 525 7 1 ?f! » J4 
2. 0 12 34 1 F .+ 3 '♦ 

2.5 S9441E + 94 
1.0430 70-: + 04 
3.553164^+04 

4. 069450'* + 04 

4.5 55 Jo IE + 04 
5. 103U16t +04 
5.674493F+04 

6. 294 H6C ♦ 04 

7. 15 0 > 66 E +04 
3.3R5o97E+04 
9.24977 7E + 04 


3. 1*1 40 5* + 03 
ft.7 T .* /35c + 13 
1.33/ l2‘5t+ -)'4 
1 . tt 3 / : 64 c +04 
2. 3 Q 37321 +04 
2.76140 7E+04 
3* 2? /1 1 3E + 34 
3.573 5465 +0'4 

4. C03.9b£>04 
4.5K J .75E + 04 
4.9flOa2bE+04 
5.28*9361+04 
5. 99975 7F +84 
6. bJ9*05C+04 
7.210 1 3 2E *0 4 


1.421 j 5*"- Q 1 
1.5179-471-0 t 
1. 55 3 '? ? \ - 0 1 
1.427 H. ’5-31 
1. Z81f7'JE-9 1 
1.14742^-01 
1,02476*1-01 
9.397'.l?f-12 
8. 7*4 *7- r -02 
8. 29174VF-02 
7.83354+F-02 
7.423 39 *r- 02 
7.0 01 3401-02 

6. 6 SB *B?F-92 
b. 127722' -02 


2.4221215-31 
7.55195?' r-31 
2.42 7*. ’756-0 L 
2.l6vH4E-11 
l.88 <l 911F-9l 
1.62 3i4 3E-’)l 
1. 354i303E-11 
1.UC256E-U 
9. 65 G ^970 0 2 
3.55 5*, HE -02 
7.687324E-02 
6. 814 76 IE -02 
5. 5123145-02 
4.507S32E-02 
3.293747E-02 


3. 92V682E-02 
1.052265E-01 
8.O85124E-02 
7.270322E-G2 
5.06748 7E-Q2 
4. 767147E-02 
3.271347E-02 
1.70 I249E-02 
8 .6&1205E-03 
2.586667E-03 
-1.462241E-03 
-6.06131 8E-03 
-1.4 8 9525E - 02 
-2.t3t086E-02 
- 2 . 82 8975E-02 


2 .8ZD369E+02 
4.3 79277 i.+62 
6 . 664779E+02 
6. 5 2831 3 E* 02 
5 . 116689E+02 
4 .651662E+02 
4 .277078E+02 
3.935920E+02 
3. B l 35 2 SE* 02 
3. 7 8 5461 E *02 
3. 708177E+02 
3 . 646032E+02 
3 . B16B0QE+82 
7.9 17020E+02 
3.752955E+02 


4.199117E+02 

5. 144674E+ 02 
4. R47V34F ♦ 02 

6. 390 9772 + 02 
5. 724431E+ 02 
4 * 902162E * 02 
4. I1794CDE+ 02 
J.692249E+D2 
3.2957112 + 02 
3.017 811E+82 
2. 477642E+02 
2.172179E+02 
3 . 328602E ♦ 02 
2. 4959942 *02 
1.8685862+82 


1. 37 89 48 r «-0 2 
1.4653 97E+02 
l.l82655 r +»2 
3.6266364*01 
6.077473E+01 
2.505002£*tt 
-1.976776E+0I 
-2.436 7 14E +01 
-5. 178 Q64E+01 
-7.676505fc+3l 
-1.230 535F+02 
-1.473 B54E+02 
-4.881984E+01 
-1.421 026E+02 
-1 ■ 804369c *02 
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